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A  tv/o-line  radiative  transfer  model  of  a  pyrotechnic  illuminating 
flare  flame  was  formulated  and  validated.  The  model  is  capable  of 
predicting  the  spectral  radiant  flux  of  different  illuminating  flares 
from  knc-ai  system  variables  such  as  formula,  si?e,  and  ambient  pressure, 
these  having  .been  varied  ever  a  wide  range.  This  was  done  without 
introducing  assumptions  which  require  ad  hoc  modifications  of  the 
modal  to  describe  different  flares. 

Relative  radiant  power  spectra  are  presented,  of  the  flames  from 
three  different  pyrotechnic  flare  formulas  burning  at  ambient  pressures 
of  760,  6?3,  300,  225,  150,  75,  30  and  6  torr.  The  sodium  concentration 
in  the  flare  formulas  varies  by  a  factor  of  10  between  each  formula. 

The  experimental  spectra  of  these  illuminating  flare  flames  show  the 
magnitude  of  sodium  D  line  broadening  as  a  function  of  ambient  pressure 
and  sodium  atom  density  in  the  flame. 

A  set  of  theoretical  spectra,  computed  using  the  two-line  radiative 
transfer  model,  are  presented  for  comparison  with  the  experimental 
spectra.  The  correlation  between  theoretical  and  experiment?, I  spectra 
snows  that  an  LTE  radiative  transfer  model  is  useful  for  prediction 
of  radiant  power  spectra  of  magnesium-' alkali  nitrate  flares,  or, 
alternatively,  these  flares  are  a  predictable  laboratory  model  radiative 
transfer  s vs ten. 
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A  two-line  radiative  transfer  model  of  a  pyrotechnic  illuminating 
flare  flame  was  formulated  and  validated.  The  model  is  capable  of 
predicting  the  spectral  radiant  flux  of  different  illuminating  flares 
from  known  system  variables  such  as  formula,  size,  and  ambient  pressure, 
these  having  been  varied  over  a  wide  range.  Th;s  was  done  without 
introducing  assumptions  which  require  ad  hot  mod :f teat  Sons  of  the 
model  to  describe  different  flares. 

To  solve  the  transfer  equation  for  observed  radiant  intensity, 
the  flame  is  represented  by  a  model  whose  main  characteristics  are 
(a)  the  flame  is  a  homogeneous  gaseous  atmosphere  with  plane-parallel 
stratification,  (b)  the  gas  consists  of  inert  molecules  plus  sodium 
atoms  which  can  be  excited  to  the  ZP  or  '  P^  level,  (c)  there  is 
local  thermodynamic  equilibrium  governed  by  the  local  temperature, 

(d)  the  temperature  gradient  can  be  represented  by  a  parabola  whose 
vertex  is  at  the  canter  of  the  flame,  (e )  th?  dispersion  profile  and 
number  density  of  sodium  atoms  have  average  values,  inside  th,=  flame, 
that  are  independent  of  depth,  and  (f)  the  individual  line  dispersion 
profile  is  replaced  with  a  two-iine  function  to  simultaneously 
describe  the  spectral  distribution  of  both  of  the  sodium  D  lines. 

The  parameters  of  the  radiative  transfer  theory  were  supplied 
from  calculated  thermodynamic  properties  of  the  flare.  Optical 
thickness  as  s  function  of  position  in  the  flame  was  determined  using 


computed  sodium  atom  densities  and  physical  flame  size  was  obtained 
photographically.  A  flame  temperature  gradient  was  constructed 
numerically  as  a  function  temperature  in  the  flame  using  the 
computed  adiabatic  temperature  at  the  flame  center  and  the  boundary. 
The  two-line  dispersion  profile  was  constructed  as  a  function  of 
line  broadening.  The  magnitude  o*  the  broadening  was  computed 
-x  priori. 

Relative  radiant  power  spectra  are  presented  of  the  flames 
from  three  different  pyrotechnic  flare  formulas  burning  at  ambient 
pressures  of  760*  630,  300,  225,  150,  75,  30  and  6  torr.  The 
sodium  concentration  in  the  flare  formulas  varies  by  a  factor  of  10 
between  each  formula.  The  experimental  spectra  cf  these  illuminating 
flare  flames  show  the  magnitude  of  sodium  D  line  broadening  as  a 
function  of  ambient  pressure  and  sodium  atom  density  in  the  flame. 

A  set  of  theoretical  spectra,  computed  using  the  two-line 
radiative  transfer  model,  are  presented  for  comparison  with  the 
experimental  spectra.  The  correlation  between  theoretical  and 
experimental  spectra  shows  that  an  LTE  radiative  transfer  model  is 
useful  for  prediction  of  radiant  power  spectra  of  magnesium-alkali 
nitrate  flares,  or,  alternatively,  these  flares  are  a  predictable 
laboratory  model  radiative  transfer  system. 
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INTRODUCTION 
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The  purpose  of  this  research  is  to  characterize  the  mechanism 
responsible  for  the  large  luminous  efficacy  of  the  magnesium-sodium 
nitrate  pyrotechnic  flare  and  to  determine  whether  rec’dtrve  transfer 
theory  predicts  the  spectral  radiant  power  of  different  illuminating 
f,  *  with  large  variations  of  system  variables  such  as  formula,  flare 
diameter,  and  ambient  pressure. 

Preliminary  research1  showed  that  the  major  component  of  the 

emission  magnesium-sodium  nitrate  flares  is  a  continuous 

spectrum,  called  the  sodium  resonance-line  continuum.  The  emission 

occurs  in  a  broad  but  variable  region  on  either  side  of  the  sodium 

resonance  lines.  The  resonance-line  continuum  from  large  flares  at 

atmospheric  pressure  extends  for  several  hundred  angstroms  mere  or 

less  symmetrically  about  the  sodium  D  lines,  wnich  are  strongly 

reversed.  The  spectrum  can  be  characterized  by  a  parameter, 

aW  =  |a„„  -  a  j  ,  the  difference  between  the  sodium  D  line  wave- 

length  ,\p  and  the  wavelength  of  maximum  spectral  flux  density 

Xmax  iocated  on  t}ie  short  wavelength  side  of  the  sodium  Ds  line. 

The  value  of  aW  increases  with  increasing  size  of  flare  and 
a 

increasing  ambient  pressure.  The  resonance-line  continuum  is  many 
orders  of  magnitude  greater  than  car,  be  attributed  tc  any  simpfe 
consideration  of  Doppler  or  Lcrentz  effects,  suggesting  that  some 
other  mechanism,  such  as  that  provided  by  radiative  transfer  theory, 
Crt-cst  give  rise  to  the  resonance-; ine  continuum.  Understanding  the 
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ori gin  of  the  resonance-line  continuum  appears  at  this  point  to  be  the 
key  to  a  more  quantitative  understanding  of  the  radiant  flux  from 
magnesium-sodium  nitrate  flares. 

Further  research^  compared  the  energy  radiated  by  the  sodium 
resonance-1 ine  continuum  with  the  energy  of  the  flare  reaction  to 
further  characterize  the  mechanism  that  gives  rise  to  tnis  continuum. 
Additionally,  an  observed  high-re^oluvi  m  radiant  power  emission 
spectrum  from  a  *0.3  cm  diameter  flare  at  at^spneric  pressure  was 
compared  to  a  spectrum  predicted  by  radiative  transfer  of  emission 
from  sodium- 0  lines,  fhe  Hunsner  and  Kybicki  formalism'  was  used  to 
solve  the  radiative  transfer  equation.  The  fig*' 2  flame  w.-v$  represent¬ 
ed  by  a  model  in  which  the  frame  is  an  isothermal  atmosphere  and  the 
Planck  function,  Voigt  function,  and  deexcitation  probability  y 
have  average  values,  inside  the  fldtite,  that  are  independent  of 
optical  depth  ~  where  1  is  the  probability  per  collision  that 
an  emitter  will  ce  deexcitsd  by  the  collision.  The  doublet  sodium  r- 
ilnes  et  589.0  and  589.8  nm  were  taken  to  be  a  single  line  with  an 
oscillator  strength  of  unity.  This  model  yielded  reasonable  agree¬ 
ment  between  the  computed  spectrum  and  tne  experimental  spectrum  cf 
a  10.8  cm  diameter  flare  vr:th  high  sodium  atom  density  burning  at 
atmospheric  pressure.  However,  the  model  is  incapacie  of  predicting 
the  spectral  distribution  of  the  P  lines  as  a  r-esolved  doublet  with 
each  of  the  D.  and  lines  strongly  reversed  at-  line  center  as 
occurs  in  low  pressure-low  sodium  atom  numt^r  density  flames. 


a^t 


MM 


furthermore,  the  defects  caused  by  this  and  the  isothermal  assumption 
we:e  somewhat  arbitrarily  mitigated  by  treating  >  as  an  adjustable 
parameter.  To  overcome  these  deficiencies,  a  more  detailed  model  was 
needed  which  included  provision  for  a  temperature  gradient  distributed 
radially  through  the  flare  flame,  for  treatment  of  the  0  lines  as  a 
doublet,  and  for  determining  all  parameters  of  the  theory  from 
ambient  ’lace  conditions  with  no  parameters  to  be  adjusted. 

The  purpose  of  the  present  research  is  to  formulate  and  validate 
a  two-line  model  of  the  pyrotechnic  illuminating  flare  flame  capable 
of  predicting  the  spectral  radiant  flux  of  different  illuminating 
flares  from  known  system  variables  such  as  formula,  size,  and  ambient 
pressure,  these  having  been  varied  over  a  wide  range,  and  to  do  this 
without  introducing  assumptions  which  require  ad  hoo  modifications  of 
the  ire-del  to  describe  different  flares.  To  do  this  we  shall  first 
present  experimental  radiant  power  spectra  of  three  different  flare 
formulas  burned  at  eight  different  ambient  pressures,  then  describe 
the  determination  of  thermodynamic  parameters:  namely,  sodium  atom 
number  density  and  adiabatic  flame  temper  crura,  and  finally  compare 
the  experimental  spectra  with  radiant  power  spectra  obtained  theo¬ 
retically  using  a  two-line  radiative  transfer  model  of  the  flame  at 
local  thermodynamic  equilibrium  containing  a  temperature  profile 
distributed  radially  through  the  flame. 
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The  spectral  radiant  intensity  of  illuminating  flames  fluctuates 
comioerably  due  to  normal  variables  of  the  flare  manufacturing 
orocess.  Flare  combustion  also  becomes  increasingly  more  irregular 
as  the  ambient  pressure  is  reduced.  The  sodium  concentration  In 
the  flare  formula  and  the  ambient  pressure  were  chosen  to  cover  a 
wide  range  to  make  variability  of  the  flare  radiative  output  small, 
compared  to  changes  due  solely  to  sodium  concentration  and  ambient 
pressure.  The  effect  of  sodium  concentration  and  ambient  pressure 
changes  is  therefore  completely  unambiguous  even  in  the  presence 
of  relatively  large  flare  output  fluctuations. 


Measurement  Parameters 

Each  uf  three  different  illuminating  composition  formulas  was 
tested  at  8  levels  of  pressure;  namely  760,  630,  300,  225,  150,  75, 

30  and  6  torr.  For  each  pressure-formula  combination,  the  burnirig 
time,  flame  size,  a:.d  relative  spectral  radiant  power  distribution 
in  the  visible  region  were  recorded. 

The  test  flares  were  composed  of  50  g  of  a  magnesium-sodium 
nitrate-binder  mixture  compressed  into  3.3  cm  i.d,  by  5.5  cm  long 
paper  tubes,  having  formulas  shown  in  Table  K  Formula  groups  1, 

2,  and  3  are  nearly  stoichiometric  mixtures ,  the  sodium  nitrate  ir, 
groups  2  and  3  being  ,!  and  .01  of  group  1  respectively.  Stoichiometry 


•5- 


was  maintained  in  groups  2  and  3  by  addition  of  potassium  nitrate 
chosen  because  it  reacts  with  magnesium  at  about  the  same  rate  as 
sodium  nitrate  end  because  of  its  low-femissivity  in  the  neighborhood 
of  the  sodium  D  lines,  the  region  of  interest  for  these  studies. 


bepen mental  Apparatus 

Fig.  1  shows  the  experimental  arrangement.  The  test  flares  F 
were  positioned  inside  a  6  m“  vacuum  chamber,  centered  between 
windows  and  1£,  ,  the  flares  burning  cigarette  fashion  with 
the  flame  projecting  upward.  Simultaneously,  as  each  flare  burned, 
the  camera  viewed  the  flame  through  window  and  plate  glass  G 
to  record  flame  'size.  The  grating  spsctrcyraoh  was  used  to  record 
the  spectral  distribution  of  the  flame  through  aperture  a  .. 
window  B1  ,  and  mirror  .  The  burning  duration  At  was 
measured  with  a  stopwatch,  fable  I-  contains  flare  burning  time 
averages  <jt.  A  He-Ne  CW  laser  and  Ar  ion  pulse  laser  were  used 
to  maintain  alignment  of  the  spectrograph  with  the  flare  and  to 
provide  wavelength  calibration  points. 

A  3m  Jarrell  Ash  Mcdel  JA-78  spectrograph  using  a  30  micron 
entrance  slit  and  fitted  with  a  racking  camera  was  positioned  to 
view  a  3  cm  wide  by  5  cm  high  region  of  the  flame  defined  by  A 
on  the  flare  axis  centered  on  a  point  about  3  cm  above  the  burning 
surface.  The  optical  path  was  changed  90°  by  a  retractable  front 


surface  Diane  mirror  H  between  the  slits  S  and  window  W  . 

* 

Kodak  Linayraph  She! i burst  35  mm  f*;1in  with  typical  usable  range  of 
400  to  700  r.m  was  exposed  to  the  flare  for  a  known  time  period  chosen 
to  provide  film  transmittance  in  the  range  0.2  to  0.8  in  the  vicinity 
of  the  sodium  D  lines,  the  region  of  maximum  interest. 


Data  Collection 
Physical  Flams  Depth 

Each  flame  was  photographed  with  a  35  mm  camera  using  Kooak 
Plus-X  film  through  an  ND2  neutral  density  filter  at  f/4  lens  stop 
for  1/125  sec.  The  exposure  settings  and  150  :m  object  distance 
were  constant  for  all  tests.  A  grid  of  reference  marks  of  known 
spacing,  photographed  while  in  the  flare  location,  was  used  to 
establish  a  linear  scale  for  measuring  the  flame  size  recorded  on  the 
film.  The  total  physical  flame  depth  z‘  was  taken  to  be  the 
distance  bet v/een  equal  film  density  regions  at  the  flame  edge 
perpendicular  to  the  flame  axis  and  through  a  point  5  cm  from  the 
flare  surface.  The  same  film  density  was  used  in  examining  all 
photographs. 

Wavelength  Calibration 

To  calibrate  the  spectrograph  film  for  wavelength,  an  argon 


ion  and  helium-neon  laser  were  exposed  to  the  film  providing  lines 
at  476.5,  488.0,  4S6.5,  514.5  and  632.8  nm.  In  addition,  Na  D 


lines  at  588.9?  and  589.59  nm,  Na  doublet  at  568.3  and  568.8  nm, 

Na  doublet  at  615.4  and  616.1  nm,  K  lines  at  578.2,  580.2.  .'81.2, 
and  583.2  nm.  and  the  Ba  line  at  553.6  nm,  appearing  in  flare 
spectra,  were  used  as  calibration  points.  The  Ba  line  appears  in  the 
flare  spectra  as  an  impurity  originating  from  residue  of  the  composition 
used  to  ignite  the  flare.  The  ignition  composition  was  10%  boron 
and  90%  barium  chromate. 

Radiant  Power  Determination 

To  determine  the  relative  spectral  radiant  power  of  the  flare, 
it  is  necessary  to  apply  corrections  which  represent  (a)  the 
relationship  between  the  irradiance  working  standard  used  during 
the  experiment  and  an  irradiance  standard  traceable  to  NBS, 

(b)  the  spectral  characteristics  of  the  window-mirror  arrangement, 
and  (c)  the  relative  power  of  the  irradiance  working  standard  and 
tbs  flare.  These  three  corrections  appear  as  time  ratios  in  the 
expression  for  the  relative  spectral  radiant  power  of  the  flare 
$  1  at  wavelength  ), 


A 


k  E^(ti/t,°)(t//tE)(t37t,) 


(1) 


where  k  is  a  proportional ity  constant  and  Ex°  is  the  working 
standard  i madia  ice.  Each  of  the  time  ratios  is,  in  effect,  a 
calibration  factor.  The  two  times  in  a  given  ratio  are  those 
required  to  expose  the  film  to  the  same  density  for  each  of  the  two 
sources  or  source  arrangements  being  compared.  The  ratios  are 


measured  at  film  positions  corresponding  t<>  each  wavelength.  In 
each  case,  t  °  represent;  the  time  for  direct  exposure  of  the 
spectrograph  film  to  the  working  irradiance  standard,  where  xrl,2,  or 

Tne  ratio  (tyt,0)  compares  system  characteristics.  The  flare 
flux  passes  through  window  W1  and  subsequently  is  reflected  by 
mirror  M  onto  spectrograph  silt  S  as  shown  in  Fig.  1.  The 
spectrograph  film  was  exposed  for  time  t,  to  the  spectral  irradiance 
working  standard  from  position  F  with  mirror  M  in  place.  This  is 
called  system  exposure.  The  film  was  also  exposed  with  the  working 
standard  in  position  L  and  with  the  mirror  M  retracted.  This  ’s 
called  direct  exposure.  The  optical  path  length  was  254  cm  in  each  ca 
The  transmittance  of  the  window  and  reflectance  of  the 
mirror  are  taken  into  account  by  t(  .  The  comparison  of  t  to 
the  time  for  direct  exposure  of  the  film  to  the  working 

standard  is  in  effect  a  correction  for  losses  due  to  the  window  and 
mirror. 

The  ratio  (t  #/tJ  compares  the  working  standard  irradiance  to 
that  of  the  flare.  In  each  case,  the  film  was  exposed  to  the  flare 
for  a  known  time  period  t  chosen  to  provide  film  transmittance 
between  0.2  and  0.8  in  the  vicinity  of  the  sodium  D  lines.  The 
ratio  {t3°/t^}  compares  the  working  standard  direct  exposure  to 
NBS  irradiance  standard  direct  exposures. 


To  obtain  exposures  at  tne  various  times  needed  to  evaluate  each 
of  the  three  ratios,  exposures  were  placed  on  each  film  in  a  similar 
pattern.  ?ach  f , 1m  contained  one  exposure  of  the  source  or  source 
arranger, ent  to  be  evaluated*  and  multiple  time  exposures  of  a  ref¬ 
erence  source.  The  > ei ?v'ence  exposure  times  were  chosen  to  include 
tre  upper  and  lower  limit,  v*'  fits  density  for  the  source  exposure. 

To  determine  the  ratio  (t  ft  ere  value  of  i  c’  for  a 
working  standard  direct  exposure  having  'he  came  film  density  as  th*. 
working  standard  system  exposure  is  interpolated  from  the  various 
reference  exposures  of  * To  determine  the  ratio  tv/7t  ,  the 
value  of  a  working  standard  direct  exposure  time  i, c  for  a  reference 
exposure  having  the  same  film  density  as  the  flo'e  exposure  t„  is 
interpolated  from  the  various  reference  exposures  t.c.  lo  determine 
t,°/t  ,  the  value  of  t  for  an  N8S  standard  direct  exposure  having 

W  O  O 

the  same  film  density  as  the  working  standard  direct  exposure  t.,° 
is  interpolated  from  the  various  exposures  of  t^. 

To  perform  the  interpolations,  the  transmittance  r  of  each  of 
i.  exoosures  on  the  film  was  measured  with  a  scanning  densitometer. 
Transmittance  values  of  the  n  reference  source  exposures  were  used 
to  construct  a  calibration  cur'e  for  the  film  for  each  wavelength 

o  o 

interval  ax(.lA  5  ax  2  1A).  The  interval  was  determined  by  the 
resolution  required.  The  calibration  curve,  a  plot  of  film  transmittance 
against  logici  ,  where  is  the  exposure  time  for  reference 

source  n,  is  rougnly  linear  over  the  useful  range.  By  interpolation. 
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a  time  (t  ° ,  t  °  or  t,  )  is  found  wnere  tne  film  densities  of  the 

*  *-  3 

source  and  reference  exposures  are  equal.  That  time  correspond:,  to 
ti  ,  or  t°  of  the  source  exposure  respectively.  The  times  t  ° 
t  °»  and  found  in  this  manner  are  used  tc-  evaluate  Eq.  (1). 


Data  Averagino 

The  film  transmittance  data  from  the  densitometer  fluctuated  over 
a  range  of  about  10%  of  the  mean  value  due  to  film  inhomcgeneities  such 
as  grammess  and  emulsion  blemishes.  To  minimize  these  fluctuations, 
the  data  were  smoothed  by  averaging  adjacent  wavelength  positions.  In 
effect,  the  procedure  was  to  record  spectra  at  higher  resolution  than 
was  actually  needed,  then  apply  to  the  resultant  data  a  mathematical 
slip  fyticpion  which  was  wider  than  the  physical  slit  width.  In  this 
running  average  method,  the  i-th  value  of  the  transmittance  is 


f  j*i+k 


-Oi  =  !  )  t(j)|/{21c+1  )  , 

I  - 

L3S1-K  J 


(2) 


where  i  ~  k+! .  k+2,  1^3,  ,  .  .  ,  n-k  and  n  is  the  number  of  data 
points  in  toe  spectrum.  The  slit  function  parameter,  2k+l ,  is  the 
number  or  data  points  over  which  the  average  is  taken.  This  is  made 
as  large  as  needed  to  achieve  de$;red  smoothing  without  distorting 
the  spectrum. 


Data  averaging  was  kept  to  a  minimum.  Smoothing  was  performed 
on  the  spectra!  correction  data  resulting  from  the  product  of  E^°, 

(t  / t  °)  and  (t  °/t  )  using  k=5.  In  this  case  the  average  value  of  a 

!  1  3  3  J 

o 

point  was  influenced  only  by  data  within  ±5  A.  Individual  trans¬ 
mittance  curves  used  to  determine  the  ab*'-,e  ratios  were  not  smoothed. 
Transmittance  values  of  the  multiple  reference  expo.sures  of  the 

working  standard  used  to  obtain  t  °  were  averaged  using  k=4 

<. 

O 

causing  only  data  within  .*4  A  to  influence  the  averaged  value.  No 
•data  averaging  was  applied  to  any  of  the-  flare  exposures. 


Resul ts 

Relative  radianr  power  spectre,  of  typical  flares  for  each 

pressure-formula  combination  are  plotted  in  Figs.  2  and  3.  Relative 
power  spectra  of  all  the  flares  tested  during  the  present  research  are 
plotted  in  Apperju  A.  The  solid  curves  in  Figs.  2  and  3  are  the  experi¬ 
mental  data.  These  spectra  were  normalized  so  that  the  peak  value  is 
unity  for  convenience  ir.  the  first  step  of  the  theoretical  comparison. 
Spectra  were  not  obtained  for  formula  groups  2  and  3  at  6  torr  because 
the  flakes  did  not  sustain  combustion  at  this  pressure.  Group  1  flares 
at  6  torr  barely  burned.  Combustion  difficulty  was  visually  observable 
for  all  flares  tested  at  75  torr  or  less.  The  lengthening  of  the  flare 
burning  time  (decreasing  burning  rate)  with  pressure  reduction  is  shown 
in  Table  II. 


The  relative  radiant  power  of  each  flare  was  obtained  by 
numerical  integration  of  the  flare  relative  radiant  power  spectrum 
i>A'  (before  normalisation)  over  the  wavelength  interval  of  interest. 
Because  the  radiant  intensity  of  these  flames  fluctuates  considerably 
even  under  normal  conditions,  flare  radiant  power  values  are  most 
difficult  to  pin  down,  particularly  at  low  p -essures  where  combustion 
is  especially  irregular.  Nevertheless,  flare  radiant  power  values, 
relative  over  the  whole  family  of  flares,  are  plotted  in  Fig.  4  for 
comparison  with  theoretically  predicted  power  values. 

Parameter  aW  1  ,  the  difference  between  the  sodium  D 

»  £ 

resonance  wavelength  Ar  and  the  wavelength  of  maximum  spectral 
flux  density  at  shorter  wavelength  than  x  ,  was  obtained  directly 

R 

from  the  flare  radiant  power  spectrum  «j>, 1 .  Values  of  aW  ' 

A  R 

for  each  pressure-formula  combination  are  plotted  in  Fig.  5 
where  they  are  compared  with  values  obtained  theoretically. 

The  flare  spectrum  half-width  aw,'  ,  measured  directly 
from  the  flare  power  spectrum  4>  '  ,  is  plotted  in  Fig.  6  for  each 

A 

pressure-formula  combination.  The  flare  spectrum  half-width,  like 
the  radiant  power  ,  fluctuates  considerably  during  normal  flare 
burning.  For  this  reason,  representative  half-width  values  are 
difficult  to  obtain.  Furthermore,  an  ambiguity  in  the  definition  of 
aW/  arises  at  pressures  low  enough  for  each  of  the  reversed  components 
of  the  Na  doublet  to  be  resolved.  The  nature  of  the  ambiguity  is 
resolved  in  the  discussion. 
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The  physical  flame  depth  z'  of  each  flare  was  measured  from  the 
photographic  negative  of  the  flare  flame.  Values  of  z'  range  from 
6  cm  for  formula  group  1,  760  torr  to  2.5  cm  for  formula  group  3,  30 
torr,  a  rather  narrow  range  considering  the  large  range  of  experimental 
conditions. 
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DETERMINATION  OF  THERMODYNAMIC  PARAMETERS 
To  solve  the  equation  of  radiative  transfer,  it  is  necessary 
to  know  the  flame  optical  thickness  and  flame  particle  velocities 
which  govern  broadening  half-widths.  These  can  be  calculated 
knowing  values  for  gaseous  sodium  atom  number  density  in  the  flame 
N  and  adiabatic  flame  temperature  T.  As  far  as  the  radiative 
transfer  model  is  concerned,  knowing  Tq  and  NQ  are  therefore 
necessary  and  sufficient  conditions  for  solution  of  the  transfer 
equation  for  the  model  to  oe  described. 

Values  for  these  parameters  could  be  obtained  relatively 
unambiguously.  The  equilibrium  composition  of  the  combustion  species 
(mole  fractions)  and  the  adiabatic  temperature  were  computed  using 
the  computer  program  developed  by  Gordon  and  McBride.4  The  program 
uses  a  free-energy  minimization  technique  to  determine  the  dynamic 
equilibrium  flame  properties.  The  calculation  recognizes  condensed 
as  well  as  gaseous  species.  Thermodynamic  functions  such  as  specific 
heat,  enthalpy,  and  entropy  are  calculated  as  functions  of  temper¬ 
ature  for  the  reactants  and  combustion  species,  for  solid,  liquid, 
and  gas  phases.  These  are  incorporated  in  the  program  in  the  form 
of  least  squares  coefficients,  having  been  derived  mainly  from  data 
taken  from  JANAF  Thermocnemical  Tables.5 

Using  ambient  pressure,  flare  formula,  and  enthalpy  values  of 
the  reactants  as  input  parameters,4  the  adiabatic  temperature 
TQ  and  the  equilibrium  composition  were  computed  for  each  pressure- 
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formula  combination.  It  is  estimated  that  30%  of  the  heat  of  combustion 
of  the  materia1  in  the  observation  region  is  lost  through  radiative, 
convective,  and  conductive  processes.  This  loss  bears  a  reasonable 
relation  to  the  25%  reportedr  for  flares  of  larger  size. 

The  ratio  of  gaseous  atomic  sodium  mole  fraction  to  mole 
fraction  of  all  gaseous  species  is  s  ,  the  atomic  sodium  partial 
pressure  being  the  product  of  with  ambient  t-rest/ure  »sc 
number  density  of  gaseot's  sodium  atoms  fj  in  the  ^Isme  was 
computed  by  the  ideal  gas  equation 


N0  «  (P3)/RT0  ,  (3) 

where  R  is  the  ideal  gar-  constant.  Values  cf  Tq  and 

are  provided  in  Table  II  for  each  ambient  pressure- formula  combination. 
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THEORETICAL 

The  Radiative  Transfer  Equation 
Diffusion  of  radiation  through  a  gaseous  atmosphere  produces 
behavior  that  is  qualitatively  similar  to  that  of  the  observed 
resonance-line  continuum.  When  the  gases  of  the  flame  are  transparent 
only  ’n  the  extreme  wings  of  the  line  scattering  profile,  maxima 
develop  on  either  side  of  the  line  and  a  minimum  develops  at  tho 
resonance  position,  where  the  optical  depth  is  greatest."  However, 
the  broadening  observed  in  pyrotechnic  flare  flames  can  be  extra¬ 
ordinary  in  comparison  with  that  which  is  normally  treated  by 
radiative  transfer  theory. 

A  radiative-transfer  mecnanism  was  tested  by  numerical 
integration  of  the  transfer  equal  ton6  using  parameters  in  the  range 
of  those  expected  in  flares  of  widely  different  formulas  burning 
at  ambient  pressures  ranging  from  760  to  6  torr.  The  total  radiant 
intensity  I  (t)  at  frequency  v  in  a  direction  described  by 

v  r* 

v  =  cose  and  issuing  from  a  volume  element  at  optical  depth  -  is 
given  by  the  radiative-transfer  equation 


ndl  (t}/di  «  $  [I  (t)  -  S  {,)]  , 

Ju  VCt  VJ  V 


(4) 


where  51  =  cosy  is  the  cosine  of  the  angle  or' observation  with 
aspect  tr,  the  outward  norms!  to  tfce  flime  surface.  A  detailed 
derivation  and-  form?1,  integration  of  tie  transfer  equation  is 
given  in  Appendix  8.  The  optical  depth  r  is  related  to  the 


physical  depth  2  hy  ,  -  rk  dz  ,  where  k  is  the  absorptivity 

X*  Vs 

of  the  flame.  The  normalized  spectral  profile  of  the  absorption 

coefficient  ±  is  a  function  which  takes  account  cf  the  flame 

va 

line  broadening  mechanisms.  Parameter  a  will  be  defined  later. 

The  line-source  function  S  (T)  accounts  for  increments  or  decrements 

\ 

in  the  radiant  intensity  from  a  volume  element  at  opticaT  depth  T 
due  to  emitters  and  absorbers  within  that  t'olume  element.  It  is 


defined  at  a  given  frequency  by  S 


^  /k  ,  where  --  is  the 

V  ')  V 


monochromatic  volume  emission  coefficient, 

formal  'ntegration  Qf  the  transfer  equation  yields  the  expression 


r 


r=r 


‘v-1  =  exp^t„‘T,)  +!  tsv(T)^a/ii]  exp[-(T-T^)^/Wjdi-  ,  (5) 

Jr- t 

where  t  and  T„  3re  the  optical  depth  integration  limits  from 
front  to  the  'ear  of  the  atmosphere  respectively,  and  1  and 
Ivs,  are  the  spectral  intensity  at  optical  depths  and  r_ 
respectively.  Ip  order  to  solve  the  transfer  equation,  for  the  observed 
radiant  intensity,  flams  is  represented  by  the  following  model. 

(1)  ;he  flame  is  a  homogeneous  gaseous  atmosphere  with  plane- 
earailel  stratification. 

(  )  The  gas  consists  of  inert  molecules  plus  sodium  atoms 
which  can  oe  excited  to  the  "P  or  "P  level. 

Jr 

{3}  There  is  Iccel  thermodynamic  equilibrium  tTE  governs*.?  by 
the  local  temperature. 


'&&*Z£fP&iSl  J^S^C 
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(4}  Energy  exchange  by  radiation  leads  to  radiative  equilibrium. 
(5}  The  refractive  index  of  the  medium  is:  unity. 

(6)  The  radiation  is  uupolarized  when  emitted  and  remains 
unpolarized  in  its  interactions  with  flame  species. 

(7)  The  temperature  gradient  can  be  represented  by  a  parabola 
whose  vertex  is  at  ti:e  center  of  the  flame. 


(8)  Thi  absorption  profile  <?vfl  ar.d  number  density  of  sodium 


atoms  No  have  average  values,  inside  the  fiame  that  are 
independent  of  t  . 

The  form  of  tq.  {5}  has  been  simplified  for  the  present  case. 

(a)  The  observed  flux  is  that  emerging  normal  tc  tin  surface  (k-1). 

(b)  No  flux  is  incident  on  the  rear  surface  of  the  atmosphere 
U^O).  {cl  Sji)  =  B.^T’)  for  the  LIE  case.  The  Planck  function 


is 


B.(T’)  »  { 2hv’'/c2 ) [exp ( hv/ kT ‘ ) - 1  ]“ 1  , 


(6) 


where  h  is  the  Planck  constant,  c  is  the  velocity  of  light,  k 
is  the  Boltzmann  constant,  and  T1  is  the  flame  temperature  at 
flame  optical  depth  ~  .  Under  these  conditions,  integrating 


from  the  front  surface,  where  z  and  r  are  0  ,  to  the  rear 


surface  where  x^  =  ?  ,  the  total  optical  thickness,  the  monoenromativ 
emergent  intensity  is 


I 
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Theoretical  relative  spectral  radiant  power  t  ,  proportional 
to  soectral  emergent  intensity  I  °  for  a  particular  model,  was 
found  by  numerical  integration  of  Eq.  (7)  on  a  C9C  6600  digital 
computer  using  Simpson's  rule  of  2m  intervals"  described  by 


f(x)dx  = 


|  [f (a)+f 


m-l  m 

f(b>+2  £  f(xaiH4  £ 
i=l  i=' 


V’4  E  rXi-d3 

i=l 


^Sf(l!(u),  (8) 


where  m>o  is  an  integer,  :i  =  (b-a)/2m,  and  x.  =  a+ih  for 

i-0,1,  —  ,  2m.  The  Fortran  program  is  listed  in  Appendix  C. 

Each  computed  spectrum  <*>x ,  normalized  so  its  power  maximum 

is  unity,  was  plotted  for  comparison  with  the  corresponding 

experimentally  determined  flare  spectrum  $. '  as  shown  in  Figs. 

2  and  3.  The  total  radiant  power  *  of  the  theoretical  flare 

spectrum,  plotted  in  Fig.  4  for  each  formula-pressure  combination, 

was  obtained  by  integration  of  •>  over  the  spectral  frequency 

region  of  interest,  the  latter  having  been  multiplied  by  By  (Tn) 

0 

where  ^  is  the  line  center  frequency.  Parameter  AWft ,  the 
separation  between  the  wavelength  of  maximum  flux  density  and  the 
sodium  line  wavelength,  and  the  spectrum  half-width  6W*  were 

^  p 

each  measured  directly  from  the  theoretical  power  spectrum  -;-v . 

Parameters  and  £Wi  are  plotted  in  Figs.  5  and  6  respectively 

<* 

for  comparison  with  the  corresponding  experimental  parameters 
i'w  '  and  &W:  !  . 

R  A 


wmm 
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Parameters  of  the  theory  that  must  be  supplied  from  properties 
of  the  flame  are  (a)  optical  thickness  t(z)  as  a  function  of 
position  in  the  flame,  (b)  a  flame  temperature  gradient  T'(z) 
as  a  function  of  position  in  the  flame,  and  (c)  the  scattering 
profile  a  parameter  in  $  .  Relative  radiant  power  spectra 

were  computed  and  compared  with  experimental  spectra.  It  remains 
to  show  that  combining  these  values  with  the  model  are,  in  fact, 
consistent  with  properties  of  the  flame. 


Justification  for  LTF.  Assumption 
A  source  function  which  does  not  assume  LIE  can  be  expressed 
in  terms  of  the  radiation  field  as 
f* 

Sv(‘)  *  [0-y)|  0(v';v,r)  va(v')  dv'3  +  y  By(i) 

Jo 

where  the  mean  intensity  J  is  the  simple  average  of  intensity 
over  all  solid  angles,  v  is  the  frequency  parameter  out  of  the 
v‘  frequency  set,  and  the  probability  per  collision  of  collisional 
deexcitation  of  sodium, 


Y  =  Cai/{C2i  +  k*S]  ~  exp(-hw0/kTo)j"  }  , 
relates  the  rate  of  collisional  deactivation  C  to  the 

C*  1 


Einstein  coefficient  for  stimulated  emission  A 


21 


(W) 


For  the 


8  - 1 


sodium  D  I’nes6  A  =  0.65x10s  s 
21 


Hunmer0  observes  tnat  when  \ 


has  its  maximum  value  of  unity,  the  source  function  S  (,)  becomes 

V 

the  Plar.ck  function  £  ( T * ) -  i.e.,  LIE  is  valid.  The  determination 
of  C  ^  will  be  described  later.  In  the  following,  it  will  be 
shown  that  a  value  of  y  large  enough  {nearly  unity)  to  justify 
an  LTE  assumption  exists  at  both  experimental  extremes  of  ambient 
pressure  and  sodium  atom  number  density. 

In  tne  case  where  the  ambient  pressure  P  -  760  torr.  the  adiabatic 
temperature  IQ  =  ZQ&'iC.,  and  the  sodium  atom  number  density  is 
ff  =  IxlO1*  cm-'" ,  v  is  about  0.35  when  nitrogen,  a  major  flame 

Q 

species  with  relatively  low  quenching  cross  section,  is  the  only 
quenching  species  considered.  An  even  larger  value  of  >  is 
obtained  when  other  flame  species  with  larger  quenching  cross 
sections  tnan  that  of  nitrogen  are  considered.  Under  conditions  of 
high  temperature,  ambient  pressure,  and  sodium  atom  number  density, 
tne  value  of  >  is  sufficiently  close  to  unity  to  justify  the 
LTE  assumption. 

At  the  experimental  limit  of  low  pressure,  temperature,  and 
sooium  atom  number  density  where  P  =  30  torr,  T  -  2500°K, 

!).  =  4.6xl014,  and  as  the  effective  Quencher,  «  is  about  3.43. 
fi  more  realistic  flame  species  mixture,  predominantly  of  N  ,  CO,  and 
CO,,  has  an  effective  quenching  cross  section  of  40x10  ‘  cm  .  For 
tills  mixture,  the  value  of  >  is  about  O.dl.  In  order  for  the 
LTE  assumption  to  be  valid  in  this  range  of  values  of  it  is 
necessary  for  tne  terms  of  Eq.  {9}  to  balance  in  such  a  way 
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that  the  resulting  source  function  S  (i)  equals  the  Planet  function 

V 

Bm(T').  Values  of  5^(0  were  compared  with  B  (V )  by  the  Hummer* 
Rybicki  formalism.'5  Some  typical  values  of  the  source  function 
Sv(t)  at  various  optical  thicknesses  for  y  -  0.5  and  a  -  0.01 
are  listed  in  Table  III.  Except  near  the  flame  surface  where  the 
optical  thickness  is  very  small,  the  source  function  Sv  is  approxi¬ 
mately  equal  to  the  Planck  value  making  the  LTE  assumption  acceptable 
ir,  this  experimental  limit  as  well. 

In  order  to  determine  the  value  of  y  by  Eq.  (10),  the 
rate  of  collisional  deactivation  C  is  needed.  C  can  be 

equated  to  the  number  of  quenching  collisions  per  second  Q,.,{T) 
given  by  HooymayersP  as 


c:  1 


C  (T)  «  Q  (T)  =  n-  v.  o. 


01) 


where  n.  is  the  density  of  tr°  quenching  flame  species,  v.  is 
J  J 

average  relative  velocity  of  approach  of  the  colliding  species, 
and  Qj  is  the  specific  quenching  cross  section.  When  the  sodium 
atoms  undergo  quenching  collisions  in  a  mixture  of  flame  molecules 
of  different  species,  the  "wenching  frequency  is  given  by* 


Ca,(T)  -  Qa,(T)  = 


m 

\  * 

;•  n-  V.  O 
—  •  j  J  j 
j=l 


(12) 


where  j  represents  various  species  of  flame  molecules.  The 
relationship10  which  describes  v.  is 


where  is  Avagadro's  number  and  M.  is  the  molecular  vieight 
of  sodium  or  species  j  as  appropriate. 

The  ideal  gas  law  was  used  to  estimate  the  number  density  of 
quenching  flame  species  n-  .  Values  of  the  specific  quenching 

J 

cross  section  a.  between  Na  and  quenching  species  0^,  ‘rt, 

CO,  CO  and  H  Q  are  8,  34,  21  ,  41  ,  50,  and  2.2xl0~1s  cmE  respectively.11 

The  dependence  of  o,  on  temperature  has  been  reported11  to  be  not 

much  stronger  than  o.  r  T~*. 

0 


Construction  of  2-line  Voigt  Function 
Under  conditions  of  low  pressure  and  low  sodium  a turn  density, 


•t  is  necessary  to  solve  the  transfer  equation  simultaneously 
for  both  of  toe  sodium  0  lines  in  order  to  describe  spectral 


distributions.  One  alternative  for  doing  this  is  to  replace  their 
individual  dispersion  profile  ;•  in  £q.  (7)  with  a  2-line  function 
V  '  This  approximation  is  valid  because  the  0  lines  are  strongly 


ccuoied.11'  thereby  maintaining  their  relative  strengths.  A  2-line 


Voigr.  function  profile  is  therefore  applicable  to  the  entire 


pressure  and  sodium  density  range  encountered  in  the  present  research. 

The  2-iine  Voigt  function,  V  “(v),  was  constructed  numerically 
by  generating  the  function  for  each  lice  separately  using  the 
procedure  described  by  Hunner , 1  then  summing  teem. 


where  V  /  ( • )  and 

Ci  1 

the  same  value  or 


are  single  line  Voigt  functions,  botn  with 
' ,  centered  on  the  Na  and  D  lino  center  frequency 


respective^.  Oscillator  strengt!  ;  ;*  and  of  the  Na  Di  and 
0  lines  are  0,312  and  0.624  respective! v. 11  Function  normalization 


and  relative  line  strength  were  maintained  by  multiplying  V  . (v)  by 
oscillator  strengths  whose  weighted  sum  equals  unity. 


Superposition  cf  Broadening  Mechanisms 
When  the  line  is  broadened  by  several  independent  effects,  the 
distributions  found  for  the  individual  kinds  of  broadening  must  be 
superimposed.  The  mathematics  of  this  superposition  is  complicated 
by  the  fact  that  the  functions  which  represent  broadening  fail  into 
two  classes,  dispersive  and  non- dispersive.  A  superposition  of 
dispersion  functions  results  in  a  new  dispersion  function  having 
a  ha If- width  equal  to  the  sum  of  the  half-widths.  However,  if 
heterogeneous  functions  are  superimposed,  each  point  on  one  interacts 
with  every  point  on  the  other.  Such  superposition,  referred  to  as 
convolution  or  folding,  of  a  Lorentz  with  a  Gaussian  distribution 
results  in  a  single  line  Voigt  function 

f  * 

v_(w)  =  f  (v-y)  fjy)  dy  ,  (16) 

CL  t  2 

where  f  (v)  is  the  Lorentz  type  profile  and  f^(v)  is  a  Gaussian 
distribution. 1 4/ ' 6 

The  normalized  forms  of  the  three  functions  are51'16 
f  (v)  =  *’1(1  +  v2)’1  Lorentz  (dispersive) 

f3(v)  =  exp(-v2) 


Doppler  (Gaussian) 
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Voigt  ( combined)  (19) 


where  v  is  the  frequency  measured  from  the  line  center  in  units 
of  Doppler  half-width  and  a  is  the  dispersion  parameter  defined 
below. 


Voigt  Function  a  Parameter 

The  Voigt  line  orofile  that  is  used  must  take  account  of  all 
of  the  factors  which  contribute  to  the  line  width  in  the  absence 
of  radiation  transfer.  Line  broadening  mechanisms  which  were 
considered  in  the  present  research  are  (a)  natural  broadening,  a 
consequence  of  the  Heisenberg  uncertainty  principle,  (b)  pressure 
broadening  of  both  Holtsmark  and  van  der  W'aals  types,  the  result 
of  collisions  with  like  and  unlike  neutral  species  respectively, 

(c)  broadening  due  to  quenching  collisions,  and  (d)  Doppler 
broadening  due  to  the  relative  motion  of  the  radiating  systems  and 
the  observer.9'14'1 ?  Pressure  broadening  caused  by  charged  perturbers 
(Stark  broadening)  was  neglected,  the  degree  of  sodium  ionization  in 
the  flame  being  small  as  shown  by  computation  of  Na+  concentration 
using  the  thermodynamic  computer  program  developed  by  Gordon  and 
McBride4  as  described  earlier.  The  collective  effect  of  these 
broadening  mechanisms  is  accounted  for  by  the  Voigt  function  a 
parameter  obtained  from  the  relation 
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a  =  E(aa  +  AA  +  AA  +  AA  )/aa  J(1n2) 
n  l  n  q  o 


(20) 


where  ax  ,  aa  ,  AA  ,  aa  ,  and  aa  are  the  natural,  Lorentz 

N  i_  f?  ■)  li 

(unlike  particle),  resonance  (like  particle),  quenching,  and  Doppler 
line  broadening  half-widths  respectively,  each  to  be  evaluated 
separately  below.  The  value  of  a  evaluated  by  Eq.  20  for  each 
of  the  three  illuminating  flare  formulas  is  6.8,  1.3,  and  0,65  for 
Groups  1,2,  and  3  respectively  at  760  torr  and  at  conditions  of 
the  flame  center.  The  values  selected  as  the  average  value  of  a 
in  the  flame  are  1.2,  0.4,  and  0.3,  which  are  in  about  the  same  ratio 
as  the  values  at  the  flame  center. 

Values  of  a  for  other  ambient  pressures  P  were  obtained 
by  ap  =  P/760,  since  the  dominant  broadening  half- widths  aa  and 
aAr  are  linear  with  the  Ambient  gas  pressure.18 

The  equation  for  Na  D  line  natural  broadening  half-width  is 


AA  (A)  *  10'  A  *i/2irTC  =  10  *.  A  /2rc 

N  0  0  U 


(21) 


where  ~  is  the  state  lifetime  and  A  is  the  line  center 

0 

-  ° 

wavelength  (cm).  AXw  is  1.2x10"  A  for  the  Na  line. 

The  Doppler  broadening  half-width  for  the  Na  Do  line  is  giyen 

ay 


U  =  2xlOs[(2  R  ln2)®/c]A  (T  /H)2‘  , 

y  o  o 


(22) 


where  M  is  tne  sodium  molecular  weight  (g).  aa^  -  0,0477  and 
0.0474  ft  at  ?c  =  2939  and  2305CK  respectively. 


The  col  1 isional  broadening  half-width  due  to  unlike  species  is 


given  by 


AX  (A)  =  1 08 ( 2x  "/rrcKi  n  {2i,kT  [(N  /m)+(N  /M)]}  ,  (23) 

u  O  c  0  A  A 

where  <?  is  the  optical  cross  section  (cm3),  n  is  the  number 

L  *»■ 

density  of  perturbing  species  (cm-3),  and  H  ar.d  m  are  the 
molecular  weights  (g)  of  the  emitting  and  perturbing  species 
respectively.  Rearrangement  of  Eq.  23  leads  to 

o  £ 

AX  (A)  =  10e(x  2/c)o  n  (8kT  N  Ary)  ,  (24) 

t  0  L  2  0  A 

where  u  *  [(mM)/(mfM)]  ,  the  reduced  mass.  If  the  mean  relative 
velocity  of  colliding  species  (cm  s"1)  is 


v.  =  (8kT  N  /vv)  , 

J  0  A 


where  j  is  the  species  index,  Eq.  24  can  be  written 


AX  (A)  -  1 0s ( X  2/c )q  n  v4 

L  O  '  U  2  J 


If  all  non-sodium  gaseous  species  are  taker,  to  be  perturbing  species, 
an  upper  limit  is  n  =  N  -  N  ,  where  N  (cm-3)  is  the  total 

2  3  O  3 

number  density  of  the  gaseous  species  in  the  flame  computed  from 
the  ideal  gas  law  at  T  .  At  760  torr  and  using  o  =  60x10" 16  cm8 

v  o  3  L 

(reported  by  Hofmann  and  Kohn13},  ax  -  0.023C,  0.0371,  and  0.0386 

L 

for  formula  groups  1,  2,  and  3  respectively. 


o  C 


Tne 

between 


t\\ 


resonance  broadening  half-width,  result 'no  from 
like  particles  at  high  gas  densities  is"0 


(A)  --  [(3^0*}e3  f  H 


3'5]/ (4n'!c*rae) 


interaction 


(27) 


*X  »<r 

where  e  is  the  elementary  char  re  (cm'  q"  s"1)  and  m  is  the 
electron  rest  mass  (o'.  Broadening  duo  to  resonance  interaction  of 
like  particles  is  dependent  on  the  inverse  cube  of  molecular 
separation  whereas  toe  broadening  by  two  unlike  particles  under 
«.hs  same  interaction  force  has  an  inverse  sixth  power  dependence 
on  the  molecular  separation.’^  Resonance  hruadening  per  atom  is 
therefore  very  large.  At  760  torr  and  f  =1.0  for  the  Na  doublet, 
calculated  values  of  a>  are  0.442,  0.0481,  and  0.00486  A  at 
N.  =  ^ . 01  xl 01 0  ,  1.1 0x1 01  ,  and  l.llxlO10  cnT?  for  formula  groups 

1,  2 .  and  3  respectively. 

Hooymayers J  describes  the  quenching  process  as  a  shortening  of 
the  radiative  lifetime.  The  equation  for  quenching  line  haif-width 

i  t, 


“V'A)  J  10*-ee[Aei  +  C£i(T)]/2tc  . 


(28) 


Calculated  values  of  are  0.0014.  0.0022,  and  0.CC27  A  at 

760  corn  for  formula  group?  1,2,  and  3  respectively. 


r  K3\3^  ^WW'S*' 


rA*s&*vwf*. 
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Optical  Thickness 

To  evaluate  Eq.  7,  the  optical  thickness  at  physical  depth  z 
of  the  sodium  resonance  tine  is 


t(z)  =  (hv  /(4„  &V  ))(N  B  -K  B  )z 

o  o  1  12  a  21 


(29) 


where  B  anci  B  ,  the  Einstein  absorption  and  induced  emission 

13  3* 

coefficients,  are  Total.  The  number  density  fi  and  N  of 

1  2 

sodium  atoms  in  the  lower  (WS)  state  and  ‘‘pper  (2P)  state, 
respectively,  have  average  values  inside  the  flame  independent  of 
flame  depth  2 (cm)  .  The  Doppler  half-width  iy  is  in  units  of 

D 

frequency.  N  and  N  are  related  co  N  by  the  Boltzr«inn 

f  2  O 

factor,  exp(-hv/kl  )  where  N  =  fi  and  N  is  negliqib1"  in 

O  •''I  s' 

wmnu/ison  to  N  .  Total  optical  thingness  ?  is  obtained 
from  Eq.  29  when  z  equals  the  tefe!  physical  flame  depth  z!, 
i.e.  t  =  t  ( z 1 ) . 


Radial  Temperature  Proille  in  Flame 
The  mass  flow  of  a  flare  flame  is  characteristically  jlong 
the  flame  axis.  The  optical  path  of  interest  in  this  research 
is  along  the  flame  radius,  perpendicular  to  the  mass  flow.  The 
radial  temperature  gradient  T!(z),  along  the  optica!  path,  is 
needed  for  evaluation  of  the  radiative  transfer  equation  whore 
T'(z)  appears  as  a  parameter  of  the  Planck  function. 
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There  a^p  two  types  zf  flame  boundaries  =>f  interest  here. 

The  phy ; real  boundary  is  defined  as  the  location  the  interface 
between  the  ambient  air  and  the  flame  medium  containing  sodium 
a terns.  The  optical  boundary  is  defined  as  the  location  in  the 
flame  medium  where  luminescence  ceases.  These  boundaries  are 
not  far  apart  physically.  The  temperature  of  the  physical  boundary 
must  be  below  that  of  the  optical  boundary. 

The  radiative  transfer  calculation  is  sensitive  to  the  temperature 

profile  and  physical  boundary  temperature.  Parabolic  temperature 

profiles  have  been  reported21  for  high  temperature  media  whose 

composition  can  be  likened  to  the  flare  Flame.  Lewis  and  von  Elbe22 

reported  that  sodium  D  line  emission  disappears  below  1775eK,  providing 

a  temperature  estimate  for  the  outical  boundary.  Lowie~7‘  assigned  a 

tempera;uro  of  ?£Q0°K  to  the  physical  (outer)  boundary  of  a  discharge 

in  sodium  vapor.  Based  on  ths  above  information*  an  approximately 

parabolic  temperature  gradient  T’(z')  was  constructed  numerically  to 

simulate  the  radial  temperature  gradient  in  the  flare  rlame.  A  l?00°k 

physical  boundary  temperature  was  used.  The  temperature  at  the  vertex 

of  the  parabola*  coincident  with  the  flame  center  (z'/2),  was  assigned 

equal  to  T  *  the  ad'abatic  temoerature. 
c 
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DISCUSSION 

Theoretical  ❖  >  and  experimental  $  1  relative  radiant  power 

spectra  are  plotted  in  figs.  2  and  3  for  the  three  formula  groups 

at  8  levels  of  ambient  pressure.  Visual  comparison  of  *x  with 

$x‘  for  each  pressure* formula  combination  shows  that  the  distribution 

computed  from  theory  agrees  quite  well  with  the  experimental 

distribution.  A  more  detailed  comparison  can  be  made  by  considering 

relative  values  of  three  parameters  which  serve  to  characterize 

the  spectrum.  These  are  AW  ami  nti  ’  the  distance  from  the 

sodium  D.,  line  wavelength  to  the  wavelength  of  maximum  flux  density, 

AWi  and  AWj'  the  spectrum  half-width,  and  $  and  v*  the 
2  2 

relative  radiant  power  of  the  sodium  D  line  emission,  the  superscript 
prime  denoting  the  experimental  parameter. 

Some  additional  features  appear  in  the  experimental  spectra 
shown  in  Figs.  2  and  3.  These  are  {a)  Mg  I  3p3P  -  4$’S  transitions 
at  51b. 7,  517.3  and  518.4  nm,  (b)  Nal  3p~P  -  4dsD  and  3p~P  -  bs2S 
transitions  at  568.3,  568.8,  615.4  and  616.1  nm.  (c)  KI  4p‘P  -  7s~S 
and  4pcP  -  5d2D  transitions  al  378.2,  580.2,  581.2  and  583.2  nm, 

{d )  a  diffuse  band,  about  5nm  wide,  with  maximum  near  574  nm 
tentatively  assigned2*  as  being  due  to  emission  between  an 
upper  levol  bound  state  to  a  lower  level  repulsive  state,  and 
(e)  a  Ss’S  -  6p’P  transition  at  553.6  nm,  barium  being  ou 
iinpurity  regaining  -Prom  the  borcn-barium  chromate  composition  used 
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to  Ignite  the  flares.  No  attempt  was  made  tc  predict  these  features 
theoretically.  Onij  sodium  D  line  emission  was  treated  in  the 
radiative  transfer  model. 

Careful  visual  examination  of  the  radiant  power  spectra  in 
figs.  2  and  3  reveals  that  at  line  center  and  in  the  region  between 
the  two  sodium  D  lines,  theory  predicts  less  radiant  power  than 
that  observed  experimentally,  particularly  when  No  is  greater 
than  10,;.  This  difference  is  expected  because  the  model  uses 
average  N0  for  all  flame  regions  and  does  not  take  account  of 
substantial  depletion  of  No  (hence  f^}  near  the  flame  boundary 
as  sodium  atoms  react  with  air  to  form  Na^O  and  Na_0„.  In  the 
Tine  wings,  the  experimental  power  tends  to  be  greater  than  that 
theoretically  predicted  i,i  several  cases  particularly  at  higher 
N()  because  of  continuum  radiation  from  condensed  flame  species 
such  as  solid  magnesium  oxide  (smoke). 

Parameters  AW,  for  the  theoretical  and  aW, 1  for  the 
*  * 

experimental  spectra  are  plotted  in  Fig.  5.  Except  for  the  group  1 
formula  at  5  torr,  where  fla^e  combustion  was  quite  irregular  making 
the  experimental  value  doubtful,  there  is  good  agreement  between 
the  two  values. 

The  theoretical  iWj.  and  experimentally  measured  AW^1 
a  a 

half-widths  of  the  radiant  power  spectra  are  plotted  in  Fig,  6  for 
comparison.  Each  of  the  continuous  solid  lines  represents  the 
theoretically  determined  half-width  of  sodium  resonance  lines 
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taken  as  a  doublet.  At  these  pressures  and  for  these  formula 
groups.,  the  individual  sodium  0  lines  overlap  to  such  a  large  degree 
that  a  half-width  of  the  single  line  cannot  be  determined 
from  the-  spectra.  The  dashed  line  and  the  individual  points  in 
Fig.  6  correspond  to  spectra  in  which  half-width  values  or  less 
than  10  angstroms  are  observed.  In  this  region,  overlap  between  the 
sodium  D  lines  is  small  enough  to  require  determination  of  the  single 
line  half-width.  Although  parameter  AVfj„'  fluctuates  considerably 
while  the  flare  is  burning  making  n  representative  value  difficult 
to  obtain,  there  is  good  agreement  in  all  cases  between  the  theoretical 
and  experimental  values  plotted  in  Fig.  6. 

The  theoretical  $  and  experimental  $'  total  radiant  power 
emitted  in  the  region  of  /ne  sodium  D  linos  are.  plotted  in  Fig.  4 
as  a  function  of  ambient  pressure,  ■?  and  both  decrease  with 
frassure,  but  the  difference  between  the  $  and  t’  increases 
as  the  pressure  decreases.  Additionally,  substantial  dispersion 
of  the  experimental  data  is  evident.  Three  factors  contribute  to 
the  observed  differences.  First,  combustion  irregularity  and  fluctu¬ 
ating  emissive  flux  were  visually  observable  for  all  flares  tested 
at  75  torr  less.  These  difficulties  became  more  apparent  as  the 
pressure  was  reduced.  Secondly,  during  the  low  pressure  experiments, 
tne  intense  radiative  zone  of  the  "lame  was  visibly  displaced 
outside  the  region  viewed  by  the  spectrograph,  resulting  in 
low  power  values  foi  <$>'.  Finally,  as  discussed  earlier,  LIE 
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approximation  is  least  acceptable  at  the  lower  pressures.  The 
increasing  deviation  from  LTE  as  the  pressure  is  reduced  makes 
assignment  of  the  Planck  value  to  the  peak  power  increasingly  unreliable, 
leading  to  overestimates  of  power  values  for  Even  in  the  presence 
of  large  ~lare  output  fluctuations  and  increasing  differences  between 
$  and  4>‘  at  the  low  pressures,  the  agreement  between  experimental 
and  theoretical  data  is  acceptable. 

In  summary,  it  has  been  shown  that  the  spectral  radiant  power 
distribution  of  a  pyrotechnic  illuminating  flare  flame  can  be 
predicted  by  a  two- line  radiative  transfer  model  which  has  been  described. 
This  can  be  done  without  introducing  assumptions  which  require  ad  hoc 
modifications  of  the  model  to  describe  different  flares.  Known  system 
variables  such  as  flare  formula,  flare  size,  and  ambient  pressure  are 
the  necessary  and  sufficient  input  needed  for  the  theoretical  prediction. 
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TABLE  I.  Flare  Formulations 


Formula  Groups 

Ingredients _ 1 _ 2 _ 3 


Magnesium 

44.0a 

40.4 

40.04 

Sodium  nitrate 

51.5 

5.15 

0.515 

Potassium  nitrate 

— 

49.95 

54.945 

Epoxy  binder  mix 

4.5 

4.5 

4  5 

c Percent  by  weight 
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Table  II  .  Flare  mean  burning  time,  computed  number  density  of 
sodium  atoms  in  flame,  and  computed  adiabatic  flame  temperature 
of  3  flare  formulas  at  3  ambient  pressures. 


Flare  Formula  Groups 

Pressure 

Quantity 

1 

2 

3 

760 

torr 

IF  sec 

28 

28 

26 

N  cm'3 

0 

l.OlxlO18 

l.lOxlO17 

1.11x10’® 

T  Kelvin 

0 

2939 

2905 

2904 

630 

torr 

IT 

31 

30 

33 

N 

0 

8.46xl017 

9.20x10’® 

9.26x10 

T 

2920 

2Gd7 

2886 

At 

30 

35 

37 

300 

torr 

N 

0 

4. 08x1  O’7 

4.44x10’® 

4.47x10- 

T 

0 

2842 

2816 

2815 

225 

torr 

It 

38 

39 

N 

0 

3.C8xl017 

3. 34x1 01 6 

2. 37x1 O’® 

I 

0 

2812 

2788 

2787 

150 

torr 

It 

35 

41 

44 

N 

O 

2.07xl0’7 

2.24x10’® 

2.26x10’* 

T 

c 

2770 

2748 

2747 

75 

torr 

It 

37 

45 

47 

N 

1.05x10”’ 

1 .14x1 016 

1  J5x10’s 

T 

I 

0 

2698 

2680 

2679 

30  ' 

torr 

It 

52 

59 

66 

N 

0 

4.29x!0’K 

4.65x1 01 5 

4. 65x10’ 4 

T 

O 

2606 

2592 

2591 

It 

75 

6 

torr 

N 

0 

8.90xl015 

9.63xl014 

9.70x10’ J 

T 

2453 

2443 

2442 

3 
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TABIE  III.  Vaiues  of  the  source  function  S(t)  as  a  function  of  depth 
■t  when  y  -0.5  and  a  =0.01. 

i 

0.7070  0.0  front  of  flame 

0.7283  0.1 

0,7919  0.5 

0.8098  1.0 

0.9540  5.0 

0.9802  l.QxlO1 

0.9951  5.0x10’ 

0.9373  1.0x10s 

0.9987  5.0x10s 

0.9990  1  0x10s 

0.9994  5.0x10s 

0.9996  l.OxlO4 

0.9997  S.OxlO4  flame  center 

0.7070  1.0x10s 


rear  of  flame 
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Fig.  1.  A  schematic  of  the  experimental  set  up. 

Legend:  G  -  glass  plate,  ,  W  -  windows,  F  -  flare, 

A  -  Aperture,  L  -  irradiance  standard,  M  -  retractable  mirror, 
and  S  -  slit  and  shutter. 


laser 

/ 


Fig.  2.  Illuminating  Flare  Flame  Spectra  for  formula  groups  1,  2,  and 
3  at  4  levels  of  ambient  pressure.  Theoretical  relative  radiant  power 
values  are  indicated  by  boxes  (*»).  Experimentally  determined 
relative  radiant  power  values  are  shown  by  the  solid  line. 
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Fig.  3.  Illuminating  Flan*  Flame  Spectra  for  formula  groups  1,  2, 
and  3  at  4  levels  of  ambient  pressure.  Theoretical  relative  radr.nt 
power  values  ^  are  indicated  by  boxes  {«*}.  Experimentally  determined 
relative  radiant  power  values  sy  are  shewn  by  the  solid  line. 
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Fig,  4,  Theoretical  <p  and  experimentally  measured  $'  flare 
relative  radisnt  power  as  a  function  of  pressure.  <{>  values  are 
shown  for  formula  groups  1,  2,  and  3  by  the  solid  lines  and 
yalues  are  indicated  by  0>A*  and  X  for  formula  groups  1,2, 
and  3  respectively. 
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Fig.  5.  Theoretical  AW  and  experimentally  measured  aW  *  widths 

P 

between  sodium  resonance  line  center  and  wavelength  of  maximum 
spectral  flux  density  as  a  function  of  pressure.  aW  values  are 

R 

shown  for  formula  groups  1,  2,  and  3  by  the  solid  lines  and  4W 
values  are  indicated  by  Q  >  A  ,  and  X  for  formula  groups  1 ,  2,  and 
3  respectively. 


solid  line,  the  single  points,  and  the  dashed  line,  aW,'  values  are 

£ 

shown  by  Q  ,  fa  ,  and  X  for  formula  groups  1,  2,  and  3  respectively. 
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APPEND  JX  A 


RELATIVE  POWER  SPECTRA  OF  ALL  FLARES  TESTED 
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Figure  H  ,  Relative  power  spectra  of  test  flare  72  .  formula 
group  i,  burned  at  760  ton*  ambient  pressure.  The  top  tv/o  spectra 
are  normalised  with  the  peak  value  equal  to  unity.  The  login  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44.0c  magnesium,  51.5%  sodium  nitrate,  and  4,5% 
binder. 
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Figure  A?  .  Relative  power  spectra  tes*  flare  75  ,  formula 
group  j  turned  at  760  torr  ambient  pressure.  The  top  two  cpectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  ns*-e  formula 
group  1  contains  44.0-  magnesium,  51.5%  sodium  nitrate,  and  4.5% 
binder. 
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Figure  A3  .  Relative  power  spectra  cf  test  fTare  145  ,  ?o*Tffj1c 
group  1,  burned  at  760  torr  ambient  pressure,  "hr  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  u^-ty.  The-  log,-,  of 
the  spectral  power  is  plotted  tw  the  bottom  spectra.  Flari  formula 
group  1  contains  44. 0^  magnesium,  51. 5%  sodium  nitrate,  and  4. 
binder. 
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Figure  A4  »  Relative  pcwer  spectra  of  test  flare  2QA  ,  formula 
group  1,  burned  at  630  tor r  ambient  pressure.  Tne  top  two  spectra 
^re  normalized  with  fcb*  peak  value  equal  to  unity.  The  log,0  of 
tu*  spectral  power  Is  plotted  in  the  bottom  spectra, 
grci'r  1  contains  44.0"  magnesium*  51  -5%  sodium  nitrate 
binder.  : 
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Figure  A5  .  Relative  power  spectra  cf  test  *lare  20B  ,  formula 
group  1>  burned  at  630  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  eguai  to  unity.  The  log,0  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44,0*  Gagnes iym,  51.5%  sodium  nitrate,  and  4.51 
cinder. 
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figure  A6  .  Relative  power  spectra  of  test  flare  20C  ,  rormula 
group  1,  burned  at  630  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  pest  value  equal  to  unity.  The  log10  of 
the  spectral  ^ower  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44. Oi  magnesium,  51.5?  sodium  nitrate,  and  4.SS 
cinder. 
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Figure  A7  .  Relative  power  spectra  of  test  flare  115  ,  formula 
group  1,  burned  at  3QQ  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44.0%  magnesium,  51.5*  sodium  nitrate,  and  4.5% 
binder. 
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Figure  AS  .  Relative  power  spectra  of  test  flare  115  ,  formula 
group  1,  burned  at  300  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log,r  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  FI  are  formula 
group  1  contains  44.0%  magnesium,  El, 5%  sodium  nitrate,  and  4.5% 
binder. 
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Figure  A9  .  Relative  power  spectra  cf  test  flare  46  ,  formula 

group  1,  burned  at  1 5C  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log,c,  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formuia 
group  1  contains  44.0%  magnesium,  51.5%  sodium  nitrate,  and  4.5% 
binder. 
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Figure  AiO.  Relative  power  spectra  of  test  flare  50  ,  formula 

group  1,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normal  izea  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  Is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44.0$  magnesium,  51.5$  sodium  nitrate,  and  4.5$ 
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Figure  All.  Relative  power  spectra  of  test  flare  54  ,  formula 

group  1,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log1c  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44.0%  magnesium,  51.5%  sodium  nitrate,  and  4.5% 
binder. 


Figure  AT 2 .  Relative  power  spectra  of  test  flare  113  ,  formula 
group  1,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  1og,c  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44.0%  magnesium.  51.5%  sodium  nitrate,  and  4.5% 
binder. 
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Figure  A13.  Relative  power  spectra  of  test  flare  128  ,  formula 
group  1,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  c  Tctra.  Flare  formula 
group  I  contains  44. OX  magnesium,  51.5%  sodium  nitrate,  and  4.5% 
binder. 
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Figure  A14.  Relative  power  spectra  of  test  flare  123  ,  formula 
group  1,  burned  at  75  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log1c  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44.0%  magnesium,  51.5%  sodium  nitrate,  and  4.5% 
binder. 
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Figure  A15.  Relative  power  spectra  of  test  flare  123  ,  formula 
group  1,  burned  at  75  to rr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44. OS  magnesium,  51.5%  sodium  nitrate,  and  4,5% 
binder. 
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Figure  A16.  Relative  poy/er  spect>'?  of  test  flare  22 A  >  formula 
group  1,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  Kitn  the  peak  value  equal  to  unity.  The  lcg10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  !  contains  44. OS  magnesium*  51.5*  sodium  nitrate,  and  4,52 
binder. 
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Figure  A17.  Relative  power  spectra  of  test  flare  228  ,  formula 
group  1,  burned  at  30  torr  ambient  pressure.  The  top  t>*0  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  1 09 , P  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44;0%  magnesium,  51.5*  sodium  nitrate,  and  4.5:t 
binder. 


Hti-A.'JVE  RADIANT  POV,£P  WELAtIVf  RADIANT  P»,W£R 


c 

0 

$ 

£• 

$ 

9 


^  **<  *  >W‘MJ  #»tf  j^s 


h 

’■ « 

A 

i\ 

I 

a  *  s 

( \ 
i\ 

* 

H 

1  \ 

:  t 

e  • 

*  1 

s 

! ; 

1 

!  J 

o  5  3 
a. 

?  ! 

J  1 
i  i 

h 

J  ! 

f  1 

1  !  1 
!  1 
i 

t  ? 

1 

*;  i  \ 

l  1 

l  1 

£;c 

|  » 
i  i 

!  !  1 
lit  \ 

' 

III  l 

5 

-1  B  s 

J  i 
!  1 

*!!  i 

5  «  n 

r 

y 

ill  1 

u 

H  l 

I  i 

i\l  s 

< 

d5  s 

/  1! 

ri  ’  j 

fj  1} 

l|  \ 

?w^V 

1  *  '^v 

QL 

:  2 

C  ! 

C  5 

\ j 

*|  • 

?  ^VyA/y 

5»  SiS  S3?  MJ  5»S  }M  5*!  5»  5!i 

WAVELENGTH  -  MJCRONS 


*65  SSC  51J  533  583  i?J  53!  532  SS3 

WAVELENGTH  -  MICRONS 


5 


;«g  ra  -*»«< 

. 


2 


\ 

l 

j 

I 


_*/■- ; 


i 


■  l 


j 

| 


***  «s*  U'  51?  3K  55'.  5*5 

WVfcLENGW  '  KI3&S 


5  r 


4 


| 


,  « 
s  * _ _ _ _ _ j 

SS  a’  35  5»  65;  S3'  5?J 

WSvElENGTw  '  MICRONS 


Figure^A'8,  Restive  power  spectra  of  test  flare  34  »  formula 

group  1,  burned  at  30  torr  ambient  pressure.  The  top  two  ~pectra 
are  normalized  with  the  peah  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  1  contains  44 .02  ^agr;eeiuml  51.5%  sodium  nitrate,  and  4.5* 


binder . 
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Figure  *20.  Relative  power  spectra  of  tesc  flare  96  ,  formula 

group  ],  burnfeu  at  30  tqrr  ambient  pressure*  The  top  two  spectra 
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Figure  A25.  Relative  power  spectra  of  test  flare  74  ,  formula 

group  2,  burned  at  760  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log;o  of 
the  spectral  powe-  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  4u.4«  magnesium,  5.15U  sodium  nitrate,  49.95S 
potassium  nitrate,  and  4.5?  binder. 
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Figure  A26.  Relative  power  spectra  of  test  flora  144  ,  formula 
group  Zy  burned  at  760  torr  amb^erft  pressure,  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log-io  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4  magnesium,  5. j5*  sodium  nitrate,  49.95? 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A27.  Relative  power  spectra  of  test  flare  144  ,  formula 
group  2,  burned  at  760  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  iog-,0  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5.15%  sodium  nitrate,  49.95% 
potassium  nitrate,  end  4.5%  binder. 
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Figure  A28.  Relative  power  spectra  of  test  flare  261A,  formula 
group  2 5  burned  at  630  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium.  5.15%  sodium  nitrate,  49 - 95% 
potassium  nitrate,  and  4,3%  binder. 
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Figure  A29.  Relative  power  spectra  of  test  flare  261B,  formula 
group  2,  burned  at  630  torr  ambient  pressure.  Th?  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log-  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  f^mula 
group  2  contains  40. 4%  magnesium,  5.15%  sodium  nitrate,  4°. 95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A30.  Relative  power  spectra  of  test  flare  262A,  formula 
group  2,  burned  at  630  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  ijnesium,  5.15?-  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.5'':  binder. 
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Figure  A31 .  Relative  power  spectra  of  test  flare  2628,  formula 
group  2,  burned  at  630  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  aqua1  to  uni  ty,  The  log,c  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5.15%  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A32.  Relative  power  spectra  cf  test  flare  262C,  formula 
group  2,  burned  at  630.  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5.1^  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A33.  Relative  power  spectra  of  test  flare  117  v  formula 
group  2,  burned  at  300  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  !cg,c.  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40. 4i  magnesium,  5.15?>  sodium  nitrate,  49 . 9SS 
potassium  nitrate,  and  4.5t  binder. 
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Figure  A34.  Relative  power  spectra  of  test  flare  1198,  formula 
group  2,  burned  at  300  torr  ambient  oressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log,c-  or 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40 . 45;  magnesium,  5.15%  sodium  nitrate,  49.95% 
potassium  ..itrate,  and  4.5%  binder. 
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Figure  A35.  Relative  power  spectra  of  test  flare  119C,  formula 
group  2,  burned  at  300  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40,4*  magnesium,  5.15?  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A26 .  Relati'-  power  spectra  of  test  flare  149A,  formula 
•group  2,  burned  sc  2d£  torr  ambient  pressure.  The  top  two  spectra 
are  normal i zee  with  the  peak  value  equal  to  unity.  The  log10  Of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  f annul  a 
group  2  contains  40.**?  magnesium,  5.15?  sodium  nitrate,  £9.95? 
potassium  nitrate,  and  4.52  binder. 
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Figure  A37.  Relative  power  soectra  of  test  flare  149B,  formula 
group  l1  burned  at  225  torr  ambient  pressure*  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  loglC,  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40,4%  magnesium,  5-15S  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.5?  binder. 
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Figure  A38.  Relative  power  spectra  of  test  flare  149C,  formula 
group  2,  burned  at  225  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40. 4*  magnesium,  5,15%  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  f39.  Relative  power  spectra  of  test  flare  47  ,  formula 

group  2,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5.15%  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A40.  Relative  power  spectra  of  test  flare  51  ,  formula 

group  2,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  tc  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5.15%  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.5%  birder. 
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Figure  A41 .  Relative  power  spectra  of  test  flare  55  ,  formula 

group  2,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log,0  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40. 4S  magn^mm,  5.154  sodium  nitrate,  49.954 
potassium  nitrate,  and  4.54  binder. 
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Figure  A42.  Relative  power  spectra  of  test  flare  112  ,  formula 
group  2,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.41  magnesium,  5.15%  sodium  nitrate,  49,95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A43.  Relative  power  spectra  of  test  flare  125  ,  formula 
group  2,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40. 4£  magnesium,  5.151  sodium  nitrate,  49,95* 
potassium  nitrate,  and  4.5*  binder. 
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Figure  A44.  Relative  power  spectra  of  test  flare  120  ,  formula 
group  2,  ourned  at  75  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  soactral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group* 2  contains  40.4%  magnesium,  5.15%  sodium  nitrate,  49,95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A45.  Relative  power  spectra  of  test  flare  23A  ,  formula 
group  2,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5.15^-  sodium  si ^9.95* 
potassium  nitrate,  and  4,5%  binder. 
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Figure  A46.  Relative  power  spectra  of  test  flare  238  ,  formula 
group  Z,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity,  ihe  iog10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40. 43  magnesium,  5.15*  sodium  nitrate,  49,95* 
potassium  nitrate,  and  4.5*  binder. 
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Figure  A47.  Relative  power  spectra  of  test  flare  35  v  formula 
group  2,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  logio 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5.15%  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A48.  Relative  power  spectra  of  test  flare  41  ,  formula 

group  2,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5.15%  sodium  nitrate,  49,95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A^g.  Relative  power  spectra  of  test  flare  96  ,  formula 

group  2,  'urned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  i,  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5. 15%  sodium  nitrate,  49.951 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A50.  Relative  power  spectra  of  test  flare  10G  ,  formula 
group  2,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  va'fue  equal  to  unity.  The  log1Q  of 
the  spectral  power  is  platted. in  the  bottom  spectra.  Flare  formula 
group  2  contains  40.4%  magnesium,  5.15%  sodium  nitrate,  43.95% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A51.  Relative  power  spectra  of  test  fiare  104-6,  formula 
group  2,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  uni ty .  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  2  contains  40. 4%  magnesium,  5.151  sodium  nitrate,  49.95% 
potassium  nitrate,  and  4.3%  binder. 
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Figure  A52.  Relative  power  spectra  of  test  flare  68  ,  formula 
group  3,  burned  at  760  torr  ambient  pressure,  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log-io  of 
the  spectral  power  is  olotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A53.  Relative  power  spectra  of  test  flare  77  ,  formula 

group  3,  burned  at  760  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A54.  Relative  power  spectra  of  test  flare  147  ,  formula 
group  3,  burned  at  76C  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 


RELATIVE  RADIANT  POWER  RELATIVE  RADIANT  POWER 


-no- 


Flare  147 


Group  26  Torr  760 


use  m  •**•*  mvBi 


so  .515  555  .Ml  .5*5  .W  -S*  .«»  .55f  .Ml  .119 

WAVELENGTH  -  HICRCNS 


1*11  ft*  *l»t  tat.tii 


■  K>  .9 It  .913  ,«!  .593  .511  .09  .535  .HI  .f* 

WAVELENGTH  -  MICRONS 


50  .511  ,9H  .Ml  MS  .515  .IK  .Ml  ,W  .511  .519  ,S5  .511  .519  .151  .MS  ,5M  590  .Ml  .555  ,515  ,515 

WAVELENGTH  *  HJCRONS  WAVELENGTH  -  MICRONS 


Figure  A55.  Relative  power  spectra  of  test  flare  147  ,  formula 
group  3,  burned  at  760  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0*515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A56.  Relative  power  spectra  of  test  flare  147  ,  formula 
group  3,  bur.ied  at  760  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.041  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A57.  Relative  power  spectra  of  test  flare  4A  ,  formula 
group  3,  burned  at  630  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  !og10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A58,  Relative  power  spectra  of  test  flare  4B  ,  formula 
group  -3,  burned  at  630  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0,515%  sodium  nitrate,  54,945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A59.  Relative  power  spectra  of  test  flare  11  &A,  formula 
group  3,  burned  at  300  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  u»  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate.  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A6G.  Relative  power  spectra  of  test  flare  118B,  formula 
group  3,  burned  at  300  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4,5%  binder. 
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Figure  A61.  Relative  power  spectra  of  test  flare  148A,  formula 
group  3,  burned  at  225  torr  ambient'  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A62.  Relative  power  spectra  of  test  flare  148B,  formula 
group  3,  burned  at  225  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4,5%  binder. 
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Figure  A63.  Relative  power  spectra  of  test  flare  48  ,  formula 

group  3,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  egual  to  unity.  The  log-so  ^f 
the  spectral  power  is  Dlotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A64.  Relative  power  spectra  of  test  flare  52  ,  formula 
group  3,  burned  at  150  torr  ambient  pressure.  The  fop  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  pov/er  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  _3  contains  40.04%  magnesium.  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A65.  Relative  power  spectra  of  test  flare  55  ,  formula 

group  3,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40,04%  magnesium*  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A66.  Relative  power  spectra  of  test  flare  111  ,  formula 
group  3,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  34.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A67.  Relative  power  spectra  of  test  flare  126  ,  formula 
group  3,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  logi0  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.315%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A68.  Relative  power  spectra  of  test  flare  129A,  formula 
group  3,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A69.  Relative  power  spectra  of  test  flare  1P.9B,  formula 
group  3,  burned  at  150  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log-jo  of 
the  spectral  power  is  plotted  in  the  bottom  spectra  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A70.  Relative  power  spectra  ef  test  flare  121  .  formula 
group  3,  burned  at  75  torr  ambient  pressure.  The  top  two  spectra 
are  normalised  with  <;he  peak  value  equal  to.  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.042  magnesium,  0.5152  sodium  nitrate,  54.9452 
potassium  nitrate,  and  4.52  binder. 
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Figure  A71.  Relative  power  spectra  of  test  flare  124A.  formula 
group  3,  burned  at  75  torr  ambient  pressure.  T he  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  !ogl0  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A72.  Relative  power  spectra  of  test  flare  124B,  formula 
group  3,  burned  at  75  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  1  og -j ^  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A73.  Relative  power  spectra  of  test  flare  24A  ,  formula 
group  3,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  !og10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A74.  Relative  power  spectra  of  test  flare  24B  ,  formula 
group  3,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.042  magnesium,  0.515%  sodium  nitrate,  54.9452 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A75 .  Relative  power  spectra  of  test  flare  36  ,  formula 

group  3,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  tne  bottom  spectra.  Flare  formula 
group’ 3  contains  40.04?  magnesium,  0.515?  sodium  nitrate,  54.345% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A76.  Relative  power  spectra  of  test  flare  42  ,  formula 

group  3,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log-io  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54,945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A77.  Relative  power  spectra  of  test  flare  44  ,  formula 

group  3,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate*  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A78.  Relative  power  spectra  of  test  flare  97  ,  formula 

group  3,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the ’peak  value  equal  to  unity.  The  log10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  birder. 
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Figure  A79.  Relative  power  spectra  of  test  flare  101  ,  formula 
group  3,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  log-0  of 
the  spectral  power  is  plotted- in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54,945% 
potassium  nitrate,  and  4.5%  binder. 
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Figure  A80.  Relative  power  spectra  of  test  flare  1 07-1 0^  formula 
group  3,  burned  at  30  torr  ambient  pressure.  The  top  two  spectra 
are  normalized  with  the  peak  value  equal  to  unity.  The  lop10  of 
the  spectral  power  is  plotted  in  the  bottom  spectra.  Flare  formula 
group  3  contains  40.04%  magnesium,  0.515%  sodium  nitrate,  54.945% 
potassium  nitrate,  and  4.5%  binder. 
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DERIVATION  AND  INTEGRATION  OF 
RADIATIVE  TRANSFER  EQUATION 

Derivation 

The  contribution  by  any  differential  volume  element  dz  •  do 
to  the  intensity  of  a  flame,  taken  to  be  a  slab  of  plane  parallel 
stratification,  is  determined  by  the  balance  of  emission  and 
absorption  of  energy  within  the  volume  element  where  dz  and  do 
are  element  thickness  and  cross-sectional  area  respectively.  By 
convention,  the  z  axis  is  taken  to  be  normal  N  to  the  slab  and 
measured  from  z  =  0  at  some  point  outside  the  slab  toward  the 
center  of  the  slab.  On  the  other  hand,  the  direction  of  flux  flow 
s  is  taken  to  be  positive  in  the  direction  of  the  observer  outside 
the  flame  with  o  being  the  angle  between  s  and  N.  Therefore, 

s  =  -z  seco  and  ds  =  -dz  sees  (B1 } 

as  shown  in  Fig,.  B1 . 

The  increment  of  intensity  lost  by  absorption  is 

-dl  a  *  k  I  ds  (B2) 

V  V  V 

where  k  is  the  linear  coefficient  of  absorption  and  I  is 

V  V 

the  specific  intensity  of  light  of  frequency  v  incident  on  the 
rear  of  the  volume  element  along  path  ds.  The  energy  emitted 
within  this  same  volume  element  is 
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dE  -  e  dv  das  dt  ds  coso  dc 

V  v 


(B3) 


where  ey  is  the  monochromatic  volume  emission  coefficient 
through  solid  angle  dw  and  time  interval  dt.  Alternatively, 
in  terms  of  its  specific  intensity,  this  energy  is 


dv  da  dt  cosO  do 


(B4) 


It  follows  from  Eqs.  (B3)  and  (B4)  that 

dI>s  a  ev  (B5) 

and  from  Eq.  (B2)  that 


<*V /is  ■  -KK 


(B6) 


Combining  Eqs.  (B5)  and  (B6),  the  emission  and  absorption  contri~ 
butions  to  the  intensity,  leads  to 


dVds  =  -KK 


+  e 


(B7) 


Substituting  Eq.  (B1 )  into  (B7)  gives  the  following  equivalent 
expressions 


dIv/(-dz  secO)  =  -kvIv  +  , 

(B8a) 

-cos©(dIv/dz)  =  -kyIv  +  ,  and 

(88b) 

v(dlv/dz)  -  kvIy  -  ev 

(B8c) 

where  u  =  cose  . 


-1 39- 


It  is  now  convenient,  mathematically,  to  combine  the  emission 
and  absorption  coefficient  to  define  a  source  function 

5  \/K  (M> 

and  to  define  the  differential  element  of  monochromatic  optica? 
depth 

dtv  =  kvdz  .  (810) 

Dividing  Eq.  (88c)  by  and. substituting  Eqs.  (B9)  and  (BIO), 
the  differential  equation  of  radiative  transfer  becomes 

u(dl  /dr  )  =  I  -  S  .  (Bll ) 

v  v  v  v 

Introducing  the  normalized  profile  of  the  absorption  coefficient 
defined  as 

r 

*V  =  VJ  k„dv  ,  (B12! 

Jo 

the  monochromatic  optical  depth  differential  element  is  described 
by 

dty  =  4>v  dr  .  (B13) 

Substituting  Eq.  (B13)  into  (811)  leads  to  another  form  of  the 
radiative  transfer  equation, 

*(dydr)  =  4,v{Iv  -  Sv)  . 


(BH) 
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Integration 

Formal  integration  of  the  radiative  transfer  section  can  be 
completed  as  follows.  Divide  Eq.  (B14)  by  u  and  rearrange  to  get 


(dly/dt)  -  (Iy^/u)  =  -J»v*v/g  . 


(B15) 


Multiply  through  by  exp{-T4>-(/y)  to  get 


[(dlv/dt)  exp(-x<}>v/u)]  -  [(Iyi>v/n)  exp(-T$y/u)]  =  -[(S^/u)  exp(-r$v/y)]. 

(BT6) 


Since  the  left  side  of  Eq.  (BiS)  equals  d[Iy  exp{-T<^/r)]/dT  5 
substitution  leads  to 


d[Iy  exp{--.*v/u)]/dT 


-C(Sv«v/n)  exp(-T4»v/u)]  . 


(B17) 


The  integral  form,  with  t  the  limit  toward  the  front  of  the 
atmosphere  and  xo  the  li  ;oward  the  rear.  Is 


d[Iv  exp(-T9v/y)] 


T~ T 


“|  [(W)  &xp(“T^v/u)jdT 

Jt=T 

o 


(818) 


Integrating  Eq.  (B18)  gives 


rwi 

PV1  exp(-T^v/|i)]  -  [Iva  exp(-Ta*v/u)]  =  -  I  [(Sv^y/u)  exp(~T9v/y)]dT 

J T=T  (BIS) 

3 


where  I  i?  intensity  at  t  and  I  is  intensity  at  t  . 

VI  1  V3  2 

Next,  by  multiplying  thrGugh  by  expfc^/y),  and  using  the  identity 
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Cexp(-xi^^/u) j[exp(+x ^^/n)]  =  exp(O)  =  unity  ,  (320) 

we  obtain 

f  T=Ti 

*'»»  "  "v2  exPC~(Tz-Ti )4»v/y]  -j  (Sv$v/m)  exp(-T<pv/p)  expvr^^/iOdT  . 

^t=t2  (B2T>) 

8y  combining  terms  under  the  integral,  changing  the  limits  and  thus 
the  sign,  Eq,  (B21)  becomes 

(  X=X 

Jvi  =  Kz  expH^-T^^/y]  +  (Sv<?v/y)  expHx-x^/uJdx  .  (B22) 

..  T=T 

1 

£q.  (822)  shows  that  at  a  given  v  and  u  ,  the  intensity  I 

emerging  from  the  atmosphere  at  any  ^  is  equal  to  the  intensity 

Kz  indent  at  x^  attenuated  by  the  atmosphere  between  t  and 

2 

1 i  p^us  1nte9ral  of  the  source  function  incrementally  attenuated 

by  the  atmosphere  between  x  and  x 

2  1  * 

The  form  of  Eq.  (822)  is  simplified  for  the  present  case  by 
considering  (a)  only  flux  emerging  normal  to  the  surface  (u=l)  and 
(b)  no  flux  is  incident  on  the  rear  surface  of  the  atmosphere  (l  =0). 
Under  these  conditions,  integrating  from  the  front  surface,  where 
2  and  ti  are  0  ,  to  the  rear  surface  where  the  total  optical 
thickness  7  =  xs,  tha  monochromatic  emergent  intensity  is 

f~-T 

V  =  *V  •  sv  exp(-x$v)dx  . 


ST=0 


(B23) 
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APPENDIX  C 

PROGRAM  LTE4  TO  SOLVE  RADIATIVE  TRANSFER  EQUATION 
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PROGRAM  LTE4  TO  SOLVE  RADIATIVE  TRANSFER  EQUATION 

This  is  the  main  program  used  to  solve  the  radiative  transfer 
equation  (Eqn.  7)  as  described  in  the  THEORETICAL  Section.  It 
computes  the  relative  radiative  power  spectrum  of  a  pyrotechnic 
illuminating  flame  from  known  system  variables  such  as  flare 
formula j  flare  size,  and  ambient  pressure.  It  considers  only 
flux  emerging  normal  to  the  flame  surface  and  that  no  flux  is 
incident  on  the  rear  surface  of  the  flame.  It  solves  the  LTE 
case  using  the  Planck  function  as  the  source  function.  A  two-line 
Voigt  profile  as  a  function  of  frequency  is  constructed  in  the 
program  for  a  specified  value  of  the  a  parameter.  The  radial 
temperature  profile  in  the  flame  along  the  optical  axis  is  constructed 
in  the  program  as  specified  by  input  parameters. 

Only  variables  in  the  NAMELIST  statement  of  the  program  are 
needed  to  operate  the  program.  These  are: 

a.  AA  -  Voigt  parameter  a. 

b.  NFREQ  -  number  of  frequency  intervals  in  Voigt  sub¬ 
routine.  The  number  of  intervals  and  width  of  the 
interval  define  the  frequency  range  over  which  computations 
are  performed. 

c.  F  -  Doppler  half-width  in  frequency  units.  When 

F  =  zero,  F  and  OWN  are  computed.  When  F  is  provided 
as  input,  DWN  must  also  be  provided  as  input. 

d.  DWN  -  Doppler  width  in  wavenumbers. 
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e.  STEP  -  input  option  used  as  a  scalar. 

f.  STEP  *  OWN  -  distance  in  wavenumbers  between  successive 
values  of  X  in  Voigt  profile.  This  product  controls 

the  width  of  the  wavelength  interval  between  computations 

g.  OS!  -  oscillator  strength  of  D1  line. 

h.  0S2  -  oscillator  strength  of  D2  line. 

f.  OS  -  oscillator  strength  sum  of  0S1  and  0S2. 

j.  XLAM  -  wavelength  of  point  halfway  between  the  D  line 
doublet.  This  parameter  is  used  to  compute  DWM . 

k.  XLAM1  -  wavelength  of  D1  line. 

l.  XLAM2  -  wavelength  of  D2  line. 

m.  Z  -  total  physical  depth  of  flame  in  cm.  Z  must  be 
provided  in  multiple  of  .062  cm. 

n.  TEMPI  -  temperature  in  kelvin  used  to  compute  F  and 
OWN. 

o.  0EN5  -  sodium  atom  number  density  in  the  flame. 

p.  M  -  integer  in  Simpson  rule  of  2M  intervals  used  to 
perform  integration. 

q.  TORR  -  ambient  pressure  in  torr. 

r.  GROUP  -  tc.mula  identification.  Only  used  for  caption 
printing. 

s„  TEMPS  -  array  of  temperatures  of  the  flame  at  various 
depths.  A  TEMPS  value  at  the  flame  boundary  {TEMPS (1 ) ) 
and  at  the  flame  middle  (TEMPS ( { Z*8 )+T ))  must  be 


-1 46- 


provided  as  input.  Between  these  points,  TEMPS  values 
at  intervals  of  .062  cm  are  input  as  necessary  to 
construct  the  profile  desired. 

t.  PUNCH!  -  logical  variable.  If  true,  non-normal i zed 
radiant  power  spectrum  is  punched  on  cards  as  well  as 
printed. 

u.  PUNCH2  -  logical  variable.  If  true,  the  radiant 
power  spectrum,  normalized  so  that  the  maximum  equals 
PLNK2,  is  punched  on  cards.  PLNK2  is  the  Planck  value 
at  XLAM  and  TEMPS ( (2*8 )+l ). 

v.  PUNCH3  -  logical  variable,  if  true,  the  luminous 
power  spectrum,  normalized  so  that  the  maximum  is  unity, 
is  punched  on  cards. 

w.  NDUPS  -  integer  variable  which  specifies  number  of 

dupl icate  sets  of  punch  card  output  one  gets  when  PUNCH1 , 
PUNCH2  or  PUNCH3  are  true. 

x.  PL0T1  -  logical  variable.  If  true,  the  computed  power 
spectra  are  plotted  by  the  printer. 

y.  DBUG  and  DBUG2  -  logical  variables.  If  true,  intermediate 
printing  takes  place. 

A  listing  of  PROGRAM  LTE4  is  given  on  the  following  pages. 


o  o  o  c,  u  o  o  uu  ou  ouuou o  u 
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PPOGRAM  LtE4(  INPUT «  OUTPUT « TAPE5= INPUT »TAPE6~0UTPUT .PUNCH)  1 

2 

STATEMENT  NUMBEDS  USED 

NEXT  STATEMENT  NU«B£R  IS  e, 

131  '  s 

»■•»»<»)!  NOTE  »♦»*****#*  ^ 

LTE3  CONVERTED  TO  LTE4  14  SEP  72  BY  COMPUTING  LUMINOUS  POWER.  7 

SUBROUTINES  COREKT.ETEBAL  AND  BLOCK  DATA  EYES  ADOED.  8 

9 

UPDATED  12  SEP  72  TO  PUNCH  OUTPUT.  10 

11 

15  AUG  72  UPDATE  TO  NORMALIZE  TO  PLANCK  VALUE  AT  XLAN  12 

AND  CENTER  TEMPERATURE.  ALSO  GETS  INTEGRAL  VALUE  OF  SPECTRUM.  13 

18 

UPDATED  31JULY72.  PUT  IN  OSCILLATOR  STRENGTH  VA| UES  FROM  IQ 

HANS  GPIEM  BOOK  PLASMA  SPECTROSCOPY.  20 

30 

PROGRAM  LTE2  IS  A  MODIFICATION  OF  LTE1  TO  INCORPORATE  A  36 

2-LINE  VOIGT  PROFILE.  37 

PROGRAM  LTE1  DOES  SINGLE  LINE  PROBLEM  IN  THERMAL  EQUILIBRIUM.  3B 

IT  USES  ASSIGNED  (INPUT)  THERMAL  GRADIENT  (SYMMETRICAL)  TABLE.  3Q 

HE  NEED  INPUT  DATA  AS  LISTED  IN  NAMELIST  INPT1  40 

IT  IS  SET  UP  TO  DO  ONLY  SODIUM  PROBLEMS.  41 

TO  00  OTHER  LINES.  CHANGE  XMASS  IN  SUBROUTINE  DOPPLER.  4? 

43 

INTEGER  7KALF  44 

COMMON" /B6/A2 (392)  45 

DIMENSION  EINT'2000)  46 

DIMENSION  XVI 72000). XV2(2000) .XV3(2000) ,PHI1 (2000) ,PHI2(2000> .  47 

«  PHI3 (2000)  48 

DIMENSION  TEMPS (201).  XVV(2000)  49 

D I HENS I ON  JL(IOl)  50 

INTEGER  P.0  51 

INTEGER  ZHP1  52 

INTEGER  GROUP  5T 

LOGICAL  PUNCH1 .PUNCM2  54 

LOGICAL  PUNCH3  55 

LOGICAL  0BUG2  56 

INTEGER  ZP1  57 

EQUIVALENCE  (XVV.XV3) . (EINT.XV2)  58 

DATA  ZZ/IHZ/  59 

DATA  P/5/.  0/6/  60 

LOGICAL  OBUG  61 

LOGICAL  PLOTI  6? 

„  63 

NAMELIST  /INPTI/  TEHPS.M.Z.F.OS.DENS, AA.NFREO.TEMPL.DWN.XLAM  64 

°  .DBUG.STEP  65 

*  .PLOTI  66 

*  »OSI .0S2.XLAH1 .XLAH2  67 

°  .GROUP.  TORP,  D8UG2  68 

*  .PUNCH1.PUNCH2. HOURS  69 

*  .PUNCH3  70 

71 

C  PUNCH1  AND  PUNCH?  CONTROL  OUTPUT  PUNCH  LOOPS  NEAR  END  OF  7? 


uuuuuuuoouo  ouooouoooouoououuuooouoooooouou  O  o  o  ooOuooo 
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HAIN  PROGRAM.  73 

NDUPS  IS  NUMBER  OF  DUPLICATE  DECKS  OF  PUNCHED  OUTPUT.  74 

IF  PL0T1  IS  TRUE. EMERGENT  INTENSITY  IS  ( LOTTED  7S 

PRINT  CONTROL  76 

0BUG2  CONTROLS  PRINT  OF  PHI1.  PHI?.  AND  PHI3.  77 

ORUG  CONTROLS  SOME  PRINT.  78 

WHEN  DBUG  IS  TRUE.  WE  GET  A  LOT -OF  PRINT  OUTPUT,  79 

WHEN  08UC-  IS  FALSE.  WE  GET  REDUCED  PRINTING.  80 

DBUG  IS  INSIDE  THE  FREQUENCY  LOOP  AND  THUS  PRINTS  A  LOT  FOR  FACH  81 

FREQUENCY.  8? 

SETTING  OF  NFRFO  TO  A  LOW  NUMBER  WHEN  DBUG  IS  TRUE  IS  RECOMMENDED.  83 

84 

STFP*DWN  =  THE  DISTANCE  IN  WAVENUMBERS  BETWEEN  SUCCESSIVE  8S 

VALUES  OF  X  IN  VOIGT.  86 

NFREO  =  NUMBER  OF  FREQUENCY  STEPS  OF  X  IN  VOIGT  NEEDED  87 

TO  COMPUTE  EINT.  88 

NFREO  CONTROLS  DIMENSION  OF  EINT  AND  XVV.  89 

NFREO  CONTROLS  DIMENSION  OF  XVI .XV2.XV3.PHII .PHI?.  AND  PHI3.  9ft 

XVV  IS  EOUIVALENCED  WITH  XV3.  91 

XVV  IS  AN  ARRAY  FOR  ALL  FREQUENCIES  OF  THE  XV.  92 

XV3  IS  ARRAY  OF  WAVELENGTH  OF  THE  2-LINF  VOIGT  PROFILE.  93 

XLAH  =  WAVELENGTH  HALFWAY  BETWEEN  THE  TWO  LINES  94 

XLAM  IS  INPUT  IN  ANGSTROMS.  95 

XLAM1  =  WAVELENGTH  OF  01  LINE.  LONGER  WAVELENGTH  OF  THE  TWO  LINE  96 

XLAH?  IS  WAVELENGTH  OF  02  LINE.  THE  SHORTER  WAVELENGTH  OF  THE  TWO  97 

LINES.  95 

XLAM1  AND  XLAH2  ARE  8CTH  INPUT  IN  ANGSTROMS.  99 

OS  =  OSCILLATOR  STRENGTH  SUM  OF  THE  TWO  LINES.  100 

OSI  IS  OSCILLATOR  STRENGTH  OF  D1  LINE*  LONGER  WAVELENGTH  OF  THE  101 

TWO  LINES.  10? 

0S2  IS  OSCILLATOR  STRENGTH  OF  02  L;"JE.  THE  SHORTER  WAVELENGTH  103 

OF  THE  TWO  LTNES.  104 

TEMPL  =  TEMPERATURE  IN  OOPPLER  USEt  TO  COMPUTE  F  AND  OWN.  I  OS 

DENS  =  ATOM  (SPECIES)  DENSITY  IN  THE  ATMOSPHERE.  N12.  106 

F  =  DOPPLER  HALF  WIDTH  IN  FREQUENCY  UNITS.  107 

WHEN  WE  INPUT  VALUE  FOR  F.  DOPPLER  ROUTINE  IS  SKIPPED.  108 

WHEN  WE  INPUT  E.  WE  ALSO  MUST  PROVIDE  INPUT  VALUE  FOR  OWN.  109 

OWN  =  OOPPLER  WIDTH  JN  WAVENUMBERS  lift 

IF  E=0« .  OR  F  IS  OMITTED.  F  IS  COMPUTED  IN  DOPPLFR  USING  HI 

TEMPERATURE  VALUE  PUT  IN  AS  TEMPL.  112 

M  =  INTEGER  IN  SIMPSON  RULE  FOR  2M  INTERVALS.  113 

INTEGRAL  IS  DONE  WITH  SIMPSON  RULE  FOR  2M  INTERVALS  114 

AA  =  VOIGT  PARAMETER  A  115 

116 

Z  =  TOTAL  PHYSICAL  DEPTH  OF  ATMOSPHERE  IN  CM.  117 

Z  IS  INPUT  IN  STEPS  OF  .062CM.  IT  IS  EXPANDED  INTERNALLY  TO  118 

END  UP  AS  AN  EVEN  WHOLE  NUMBER.  (  er  ZA  ) 

120 

TEMPS  ARE  THE  TEMPERATURES  OF  THE  ATMOSPHERE  AT  VARIOUS  DEPTHS.  121 

THE  INOEX  OF  TEMPS  IS  THE  INTEGER  CORRf SPONDING  TO  16TJME5  THE  122 

NUMBER  OF  .062CM . STEPS  OF  7  PLUS  1.  123 

TEMPS  INOEX  =  (  7*16 J  *1  124 

EXAMPLE.  WHEN  DEPTH  =  ZERO  CM.  INDEX  =  1.  125 

WHEN  DEPTH  =  2.5  CH.  INDEX  =41.  126 

TO  MULTIPLY  BY  16 ASSURES  AN  EVEN  NUMBER  CORRESPONDING  TO  THE  127 


o  u  o  u  o  o 
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TOTAl.  PHYSICAL  DEPTH.  128 

THIS  ALLOWS  ONE  TO  INPUT  A  TEMPS  AT  INTERVALS  OF  .062  Cf  529 

TEMPS  MIDDLE  INDEX  =  <Z<*8>*1.  130 

WE  MUST  PROVIDE  TEMPS  VALUES  FOR  DEPTH  =  O.OEPTH  =  MIDDLE  (CM).  131 

AND  OTHER  VALUES  BETWEEN  0  AND  MIDDLE.  132 

TEMPS  ARE  INPUT  IN  DEGREES  KELVIN,  133 

134 

1  CONTINUE  135 

C  RESET  INPUT  VARIABLES  TO  DEFAULT  VALUE  136 

DO  2  J=1 .201  137 

TEMPS <J)=0.  138 

2  CONTINUE  139 

OWN-O «  140 

Z=4.0  141 

F=0.  142 

OENS-=0.  143 

AA=0.  144 

STEP=1 •  1*5 

DBUG=. FALaE.  146 

DBUG2=. FALSE.  147 

IHR=0  146 

HIN=0  149 

ISEC=0  150 

H=0  151 

M=200  152 

H=100  153 

H=50  154 

OS1=0 ,312  155 

OS2=0,624  156 

0S=0Si*0S2  157 

NFREQ=0  158 

NFREQ=300  159 

TEHPL=0.  160 

TEHPL=3098.  161 

XLAH=0.  162 

XLAH=5692,935  163 

XLAHl=5fl95.92  164 

XLAH2=5B89.75  165 

PL0T1=. FALSE.  166 

PL0T1=,TRUE.  167 

PUNCHl=.TRUE.  168 

PUNCH1=.FALSE.  169 

PUNCH2=, FALSE.  170 

PUNCh3=. FALSE.  171 

PUNCH2=.TRUE.  '  172 

PUNCH3=.TRUE.  173 

NDUPS=1  174 

175 

READ(P.INPTl)  176 

177 

C  PROGRAM  TERMINATES  WHEN  M=0,  173 

IF(M.EO.O)  CALL  EXIT  179 

180 

WRITE(O.ll)  181 

11  FORMAT (1H1)  182 
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CALL  DATE (DM, YE) 

C  CALL  CLOCK (IHRvHIN.ISEC) 

WRITE (Q, 63)  DM.YE, IHR.HIN, I SEC 

*  .GROUP,  TORR 

63  FORMAT  (2X*A10.A2»4X,3I3» 

*  SOX,«GROUP*I3.  F10.0,  *TORR<*> 

WRITE (G.24)H,Z.F »OS.DENS.AA»NFREQ,TEMPL»DWN»XLAN,DBUG»Sf£P 
»  tPLOTI 

24  FORMAT (/*  M  IN  THE  SIHPSON  RULE  IS*  16./ 

«  •  7.  THE  DEPTH  OF  THE  ATMOSPHERE  IN  CM  IS*  F6, 1 »/ 

*  *  F  THE  OOPPLER  HALF  WIDTH  OF  THE  LINE  IN  FREOUENCY  UNITS  IS* 

*  £12.5./ 

*  *  OS  IS  THE  OSCILLATOR  STRENGTH  SUM  OF  THE  TWO  LINES.  IT  IS* 

*  F 1 0 . 3  s  / 

»  «  DENS  THE  SPECIES  OENSITY  IN  THE  ATMOSPHERE  IN  PARTICLES/CC  IS* 

*  £12.5./ 

*  «  AA  THE  LITTLE  A  PARAMETER  IN  THE  VOIGT  PROFILE  IS*  F10.5,/ 

*  *  NFREO  IS  THE  NUMBER  OF  STEPS  OF  FREQUENCY  IN  EMERGENT  INTENSIT 
*Y.  IT  IS*  16./ 

*  «  TEMPL  THE  TEMPERATURE  IN  THE  DOPPLER  COMPUTATION  I3»F9.3»/ 

*  *  OWN  IS  OOPPLER  WIDTH  IN  WAVENUMBERS.  OWN  IS*FIG.3»/ 

*  *  XLAH  THE  WAVELEN6TH  HALFWAY  BETWEEN  THE  TWO  LINES  IN  ANGSTROMS 

*  IS* 

«  F10. 3,/  *  DBUG=*  L4,/ 

*  *  STEP  TIMES  OWN  IS  THE  DISTANCE  IN  WAVENUMBERS  BETWEEN  SUCCESSI 
*VE  VALUES  OF  X  IN  VOIGT.  STEP  =*F7.3«/ 

*  «  PL0TI=*L4  J 

WRITE (Q, 81 )0S1 .OS2.XLAHI .XLAH2 
81  FORMAT ( 

*  *  OS1  THE  OSCILLATOR  STRENGTH  OF  THE  D1  LINE  IS*  FI0.3,/ 

*  «  OS2  THE  OSCILLATOR  STRENGTH  OF  THE  D2  LINE  IS*  F10.3./ 

*  *  XLAK1  THE  WAVELENGTH  IN  ANGSTROMS  OF  THE  D1  LINE  IS*F10.3./ 

*  *  XLAH2  THE  WAVELENGTH  IN  ANGSTROMS  OF  THE  D2  LINE  IS»F10.3,/ 

*  ) 

C  INITIALIZE  VOIGT  FUNCTION  ROUTINE. 

OUMHY1-COFVOI (0UHMY2.DUMMY3J 

ZA- 16. 

Z4=Z*2A 
062=0.062 
0125=0.125 
ZP1=  IFIX (Z4) *1 
IF (.NOT, DRUG)  GO  TO  86 
WRITE (0,10)  062. (TEMPS (J> ,J=1 ,ZP1 ) 

86  CONTINUE 

10  FORMAT (*  THE  TEMPERATURES  AT  VARIOUS  DEPTHS  IN  THE  ATMOSPHERE* 

*  /«  WHERE  THE  TEMPS  INDEX  CORRESPONDS  TO  THE  PHYSICAL  DEPTH  AT  IN 
“TERVALS  OF*  F6.3*  CH.  ARE* 

«  /  (SX, 10F8.0) ) 

WRITE(0,93>  062,0125 
93  FORMAT ( 

»  *  THE  FOLLOWING  RELATES  TO  TEMPS  INDEX  */ 

*  *  INDEX* 1  MAPS  TO  OEPTH=ZERO  CH*/ 

*  *  INDEX=2  HAPS  TO  DEPTh=*  F6.3*  CH,*/ 


183 

184 

185 

186 
137 
188 

189 

190 

191 
lv2 

193 

194 

195 

196 

197 

198 
1'99 
200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 
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*  «  INDEX=3  HAPS  TO  0EPTH=*  F6.3*  CM, ETC.*  238 

*  /»  TEMPS  MIODLE  INDEX=(Z  TIMES  8>*1  *  239 

*  /)  240 

24 1 

NEXT  COMPUTE  F  AND  OWN  WHEN  NOT  PROVIDED  IN  7NPT1 .  242 

TEMPL  =  TEMPERATURE  UNPUT)  FOR  DOpPLER  SUBROUTINE.  243 

IN  DOPPLER  WHEN  2ND  ARGUMENT  =  0.  WE  FIND  F  ANO  DSN  COUTFUT).  244 

AKW  =  COMPUTED  DOPPLER  HALF  WIDTH  IN  ANGSTROMS  AT  TEMPL.  245 

F  =  OOPPLER  HALFWIDTH  IN  FREQUENCY  UNITS.  246 

OOPPLER  IS  ENTERED  WHEN  NO  VALUE  FOR  F  IS  GIVEN  IN  NAMELIST  INPT1.  247 

DOPPLER  IS  SKIPPEO  IF  WE  PROVIDE  VALUE  FOR  F.  ?4R 

OWN  =  DOPPLER  WIDTH  IN  WAVENUMBERS  249 

250 

IF (F  ,NE.  0.)  GO  TO  21  251 

CALL  DOPPLER (XL AN.  0.  TEMPL.  F.  OWN  «AHU  )  25? 

WRITE10.22)  F.  DWN.6HW  253 

22  FORMAT!//*  THE  COMPUTED  OOPPLER  HALF  WIOTH  IN  FREQUENCY  UNITS  =*  254 

*  E12.5./  255 

*  *  THE  COMPUTEO  OOPPLER  WIDTH  OWN  IN  WAVENUMBERS  =«F10.3./  256 

*  *  THE  DOPPLER  HALF  WIDTH  AHW  IN  ANGSTROMS  IS*  F1C.5./  )  257 

21  CONTINUE  258 

259 

COMPUTE  PLNKI  260 

PLNKI  =  VALUE  OF  PLANCK  FUNCTION  AT  FREQ  XLAM  AND  3098  K.  261 

CALL  PLANCK (  3098.*  XLAH.  2  .Pt.NKl  )  262 

KRITE(Q.23)  XLAH. PLNKI  263 

23  FORMAT (//*  VALUE  OF  PLANCK  FUNCTION  AT  CENTER  FREQUENCY*  FI 0.3,  264 

*  *  ANGSTROMS  IS*  E12.5./)  265 

266 

FOR  CONVENIENCE.  WE  INPUT  TEMPS  ONLY  FROM  DEPTH-0  THRU  DEPTH  267 

=  MIOOLE.  260 

WE  ASSUME  THE  TEMPERATURE  GRADIENT  IS"  SYMMETRICAL  ABOUT  THE  269 

HIDDLE.  NEXT  WE  GENERATE  THE  TEMPS  BETWEEN  MIDDLE.  ZHALF,  AND  270 

THE  REAR  BOUNDARY  OF  THE  ATMOSPHERE.  Z.  271 

272 

ZHALFc  lFIX(Z4)/2  273 

ZKPl=ZKALF.l  274 

00  13  J=l, ZHALF  275 

TEMP5 ( J*ZHP1 )  *  TEMPS(ZHPl-J)  276 

13  CONTINUE  277 

IF ( .NOT .OBUGI  GO  TO  07  278 

WRITE(Q.iO)  062. (TEHPS ( J) ♦ J=1 »ZPI )  279 

67  CONTINUE  280 

AT  THIS  POINT,  WE  HAVE  AVAILABLE  A  TEMPERATURE  GRADIENT  AT  201 

SELECTED  DEPTHS  OVER  THE  ENTIRE  ATMOSPHERE  282 

283 

HP1=H*1  284 

M2=M*2  285 

M2P1=M2*1  286 

WRI7E(Q»12)  M.H2  267 

12  FORMAT («  H  IS  THE  SIMPSON  RULE  MIDDLE  INOEX  CORRESPONDING  TO  THE  260 

*PHYSICAL  DEPTH  OF  THE  MIDDLE  OF  THE  ATMOSPHERE.  H=*I3,/  289 

*  *»  THE  SIMPSON  RULE  INDEX  CORRESPONDING  TO  THE  TOTAL  DEPTH,  2H.  OF  290 

*  THE  ATMOSPHERE  =*I4,///)  291 

WRITE(Q.IS)  292 


o  o  on  o  no 
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15  FORMAT <//»  Z  MAPS  TO  2M  AND  ZHALF  HAPS  TO  M  */  293 

«  *  ZPl  HAPS  TO  M2P1  AND  ZHPI  MAPS  TO  MPl*//>  294 

295 

FIND  THE  TEMPS  FOR  EACH  OF  THE  KISSING  DEPTH  POINTS  BETWEEN  296 

DEPTH  =  ZERO  AND  DEPTH  =  MIDDLE  (ZHALF)  297 

TEST  FOR  ZERO  TEMPS  AND  RECORD  THEIR  LOCATION  298 

K=G  299 

DO  16  J=2*ZHF1  300 

Jl=J  301 

IF (TEMPS  <J) )  50*50*16  302 

16  CONTINUE  303 

GO  TO  17  304 

50  CONTINUE  305 

K=K*1  306 

JL(K)=J1  307 

GO  TO  16  308 

17  CONTINUE  309 

NOW  ALL  THE  ZERO  TEMPS  INDEX  NUMBERS  ARE  STORED  IN  JL  310 

THERE  ARE  K  OF  THEM  BETWEEN  ZEPO  CM  AND  ZHALF  311 

312 

NEXT  WE  INTERPOLATE  BETWEEN  THE  TEMPS  AND  FIND  A  TEMPS  FOR  ALL  313 

THE  TERMS  THAT  WERE  ORIGINALLY  ZERO.  314 

IF  t  K  .£ Q.  OJ  GO  TO  35  315 

L=0  316 

DO  SI  J=2*K  317 

J1=J  318 

IF ( JL (  J)  ,EO.  ( JL { J“  D  ♦  1 > IGO  TO  53  319 

L-L*l  320 

L?1=L*1  321 

GO  TC  52  322 

53  L=L*1  323 

LP1=L»2  324 

56  DUH'O.  325 

51  CONTINUE  326 

52  CONTINUE  327 

JA  =  JL(JI-L)  -  I  328 

J8  =  JA  ♦  LP1  329 

IF (K  .EO.  1  I  JA=I  330 

IF  (K  .EO.  1  )  JB=3  331 

JBHI  *  JB-I  332 

JAPlsJA.l  333 

TOIFF  =  TEMPS (JB)  -  TEMPS (JA)  334 

TADD  a  TOIFF  /  LP1  335 

DO  54  IaJAPltJBHl  336 

TEMPS (I)  =  TEMPS (1-1)  ♦  TADD  337 

L=0-  338 

54  CONTINUE  339 

IF(K  .EO,  1  >  GO  TO  55  340 

IFCJ1  .EO.  K>  GO  TO  55  341 

GO  TO  56  342 

55  CONTZNUE  343 

IF ( .N0T.0BUG5  GO  TO  88  344 

WR ITE (0*10 >  (TEMPS (J)  « J=1 *ZP1 J  345 

88  CONTINUE  346 

C  NOW  WE  HAVE  TEMPS  FOR  ALL  POINTS  BETWEEN  Z=0  AND  ZHALF  347 


OOUO^OOOUOOU  OUfcJOOOVWOOOO  0-0  uu  oooooooo 
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34* 

349 

C  NEXT  riNO  THE  TEMPS  BETWEEN  ZHALF  AND  Z  BY  SYMMETRY  350 

DO  -18  J-  1  *  ZHALF  351 

TtMP5(J*7HPl)  =  TE.‘t°S  (ZHP1  -J)  35? 

18  CONTINUE  353 

WRITE <Q«iO)  062* (TEHPSIJ  ,J=1,ZP1>  354 

C  NOW  *£  HAVE  A  TEMPS  FOR  EACH  CM  OF  DEPTH  FROM  ZERO  THRU  2.  355 

356 

WE  INTEGRATE  BY  SIMPSON  RULE  FROH  0  TO  2M  INTERVALS.  THAT  357 

IS  FROM  INDEX  1  THRU  H2P1.  3SB 

WE  NEED  A  SUBROUTINE  WHICH  COMPUTES  A  DEPTH, ZOU  AND  A  359 

TEMPERATURE  TEMP?  FOR  EACH  INDEX  POINT  IN  THE  INTEGRATION.  360 

WE  CALL  IT  DEPTMP.  361 

IT  IS  CALLED  BY  3fc? 

CALL  DEPTMPtHfl, TEMPS, M2, Z.TEHPP.ZCH)  363 

MS  IS  THE  SIMPSON  RULE  INDEX  FOR  WHICH  WE  FIND  A  TEMP?  ANO  ZCM.  364 

TEHPP  IS  T^f 'TEMPERATURE  IN  KELVIN  AT  MB  {OUTPUT)  365 

ZCM  IS  THE  DEPTH  IN  CY  AT  MS  (OUTPUT)  366 

TEMPS  IS  THE  TEMPERATURE  INPUT  ARRAY  367 

H2  IS  2H  INTERVALS  OF  SIMPSON  RULE  36B 

369 

370 

371 

8  E  DOHDA  20  MARCH  72  37? 

OOES  SIMPSONS  RULE  FOR  2M  INTERVALS  373 

SEE  EON  6-66  OF  MOURSUNO  AND  DUBIS  374 

ELEMENTARY  THEORY  AND  APPLICATION  OF  NUMERICAL  ANALYSIS  375 

HCGRAW  HILL  1967  376 

LET  M  =  INTEGER  .GT.  ZERO  377 

LET  H  =  (B-A) /2K  378  - 

A  =  INTEGRAL  LOWER  LIMIT  379 

8  -  INTEGRAL  UPPER  LIMIT  33C 

X(I)=A*I«H  FOR  1=0,1 ...» iCM  "81 

TERM'}.  =  2«S'JM  OF  F(X(2I1)  FOR  I=i  .H-l  3S2 

TERM4  =  4°SUM  OF  F(X <21-1)1  FOR  1=1,2, . K  333 

TERH1  =  F (X { A ) >  =  F (X ( 1=0) )  384 

TKRH2  =  F(X (B) )  =  F (X ( I=2M) )  385 

INTEGRAL  =  (H/3) * (TERM1 +TERM2  *  TERH3  ♦  TERM4  )  386 

WE  NEED  TO  SOLVE  FOLLOWING  EON  AT  EACH  FREQUENCY  337 

HINT  =  PHI  *  INTEGRAL  OF  (PLANCK»EXP (~TT*PHI ) )  388 

INTEGRAL  LIMITS  ARE  A=0*  B=M2  389 

INTEGRAL  INDEX  GOES  A=l»  B=M2P1  390 

EINT  =  EMERGENT  INTENSITY  391 

PHI  =  VOIGT  PROFILE  OR  OTHER  ABSORPTION  PROFILE  392 

PLANCK  =  PLANCK  FUNCTION  AT  TEMPERATURE  CORRESPONDING  TO  DEPTH  393 

IN  THE  ATMOSPHERE  MAPPED  TO  INTEGRAL  INDEX  394 

TT  =  TOTtL  THICKNESS  AT  DEPTH  CORRESPONDING  TO  INTEGRAL  INDEX  395 

3"}£ 

C  START  tJERF  TO  GENtSATt  2-LINE  VOIGT  PROFILE.  797 

VGT=Vr>i«f  (O..AA)  393 

VGT)  . -VOIGTU0..AA)  399 

WR»TE'Q,94}  VGT • VGT 10  400 

94  FORMAT (/, 10X*  VOIGT  AT  CENTER  FREQUENCE  T£*E12.5,/«  401 

9  *  VOIGT  AT  FREQUENCY  X=10  IS*”  Ei2.5,/  5  402 


o  o 
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1 -E-8) 


69 


.32 

90 


TO 


83 

91 


73 

74 


75 

76 


76 

79 


77 

30 

72 


NFR£QH=NFRfC/2 
XWAV£NO=i ,0/<  XLAHi  « 

DO  69  J=1 .NFREQ 
X=(FLOAT { J-l-NFPEGH) )  *  STEP 
PHII (J>=( VOIGT (X.AAH  *051/05 
XVI (J5=  l.E*8  /  (XWAVENO  -  (DWN*X>) 

CONTINUE 

IFf  .NOT.  D8UG2)  GO  TO  90 
WRITE (0,32) 

FORMAT!//*  PHI 1  FOR  THE  LONGER  WAVELENGTH  LINE*) 
WRITg; (Qt47|  ( ( (XV1 1  j,  ,phii  (j) ) ,  J=1  .NFREO) ) 

CONT INUE 

XWAV£NO=1.0  /  (XLAH2  »  l.E-8) 

DO  70  J=l#NFREO 

STEP 
OS2  /OS 


X=  (FLOAT { J-I-NFREOH) )  » 
PHI2(J)= (VOIGT (X.A45 ) 


XV2 ( J) =1 «£*8  /  (XWAVENO  -  (OWN*X)> 

CONTINUE 

IF(  .NOT.  DBUG2)  GO  TO  91 
WRITE(Q,83) 

unTTrly7*  PHT?  F0R  THE  SH0RTER  WAVELENGTH  LINE*) 
WRITE(°,47) ( ( (XV2( J) ,PHI2< J) > , J=1 , NFREO) ) 

CONTINUE 

XKAV£N0=1.  /  (XLAK  *  l.E-8) 

DO  71  J=l, NFREO 

X= (FLOAT (J-l-NFREQH) )  *  STEP 

XV3 ( J) =1 .£*8  /  (XWAVENO  -  (DWN*X>) 

CONTINUE 

NEXT  AT  EACH  XV3  WE  FIND  THE  INTERPOLATED  PHII  AND  PHI? 
AND  ADD  THEM  TOGETHER  TO  GET  PHI3. 

DO  72  J=l, NFREO 
DO  73  K=1,NFREQ 
KS=K 

IF(XV2(XS)  ,GE.  XV3( J) )  GO  TO  74 

CONTINUE 

60  TO  75 

FRAC  = (A8S>  (XV3 ( J) -XV2 (XS-1 ) ) ) / (A8S (XV21K5) ~XV?(K3-1 \ > » 
PHI2  J-PHI2  (KS-1  >  *  (FRAC*;PHI2!xsTyHf2(KS-nT) 

GO  TO  76 
CONTINUE 

PHI2J=PH 12 (NFREO) 

CONTINUE 

IF (XV3 ( j)  .  (_c.  XVI  (I)  )  GO  TO  77 

DO  73  X=l, NFREO 

K5=K 

IF (XV3 ( J)  .LT.  XVI (KS))  GO  TO  79 

CONTINUE 

GO  TO  72 

K5=KS“1 

FRAC=(A8S(XV'3(J)-XV1 (KS) ) )/(AS3(XVl (KS)-XVl (KS*1 > ) j 
PHIlJsPHIl (K3)  .  «FRAC*fPKIllKS-ii-PHll(KSI») 

GO  Tn  31) 

PHIlw-PHIl  (1) 

PHI3CJ)  =  PHII J  ♦  PH 12  l 
CONTINUE 


403 

404 

405 
606 
407 

403 

404 

410 

411 
61? 

413 

414 

415 

416 

417 
41<s 

419 

420 

421 

422 

423 

424 

425 

426 

427 
423 

429 

430 

431 
43? 

433 

434 

435 

436 

437 
430 
43° 

440 

441 
44? 

443 

444 

445 

446 

447 
443 

449 

450 

451 
45? 
453 

456 

455 

456 

457 
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IF (  .NOT.  0BUG21  GO  TO  92 
WRITE (Q»84) 

84  FORMAT  (//<*  PHIS  THE  2-LINE  VOIGT  PROFILE*) 

WRITE  (9  7!  I  (  (XV3  ( J)  »PHI3t J) ) , J=l  .*;FPEO)  ) 

92  CONTXNL-" 

NOW  WE  <A  E  THE  2-LINE  VOIGT  PROFILE  IN  PHI3U). 

E4CH  I?  WE.GHTfO  FOR  OSCILLATOR  STRENGTH  AND  IS  SYMMETRICAL 
ABOUT  XL AM. 

ALL  THE  FOLLOWING  IS  INSIDE  MASTER  DO  LOOP  OVER  FREQUENCY.  X 

00  999  JJJ=1,NFPE0 
FIRST  WE  FIND  TERMl 
TERHJ  =  FUNCTION  AT  ZERO  CM  DEPTH 
TT=0 .  AT  ZERO  DEPTH 


CALL  THICK ITT.F «OS»0. »C£NS) 

PHI=PHI3< JJJ) 

EX  =  £XP(-TT«PH!> 

TEMPERATURE  =  TEMPS  ID  AT  ZERO  OEPTH 

XV  IS  ARBITRARY  FREQ  X  CONVERTED  TO  CORRESPONDING  WAVELENGTH 
IN  ANGSTROM  UNITS. 

XV=XV3!JJJ) 

TFH  =  TEMPS ( 1 ) 

CALL  PLANCK (  TEM  .Xy,  2.PLNK) 

IF ( .NOT .  DRUG )  GO  TO  25 

WRITE(Q»26)TT  *X.AA*PHI»£X«XY»XLAM.0VN»PLNK*PLNK1»XWAVEN0 

26  FORMAT!//*  TT=*fl2.5«/*  X=*£12.5/«  AA=*F10.3/«  PHI=®E1?.S/ 

*  »  EX=* 

*  E12.5/0  XV=»ri0.3/*  XLAM=»F10.3/  *  DWN=*FI0.3/*  PLNK=*E12.5 

*  /*  PLNKl=*E12.5  /»  XWAVENO  =«F10.3  ) 

WWITEI0.30) 

30  FOPMATI  *  Z=  0  CM® ) 

WRITE (Q.271 T£MPS( I ) 

27  FOFMATC*  TEMPS (1 )=*T10.3) 

25  CONTINUE 

NORMALIZE  SO  DLNK  AT  CENTER  FREQUENCY  AND  30S8K  =  1. 

PLNK  AT  XLAH  AND  30S8K  =  PLNK1 
PLNK=PLNK/PLNK) 

PLNK1  IS  COMPUTEO  EARLY  IN  THE  PROGRAM. 

TERMl  =  PtNK®EX 

IF (.NOT.  DRUG)  GO  TO  29 
WRITE (0.28)  PLNK.TERH1 

28  FORMAT l*  PLN*  FORMALIZED  =®E12.5/ 

*  30X.  *  VALUE  OF  r£RMl=*c 12.5) 

29  CONTINUE 


NEXT  WE  FIND  TFPM2 

TERM2  =  FUNCTION  AT  Z  CM  DEPTH 


458 

459 

460 
46f 
*6? 
469 
469 

465 

466 

467 

468 

469 

470 

471 
47? 

473 

474 

475 

476 

477 

478 
679 

480 

481 
48? 

483 

484 
4SS 
466 

487 

488 

489 

490 

491 
49? 

493 

494 

495 

496 

497 
495 
494 

50  0 
501 
50? 

503 

504 

505 

506 

507 

508 

509 

510 

511 

51  2 
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7.7=1 

CALL  THJCK(TT,F,0S.?7,DEN5> 

C  SAVF  TT  IN  0!J1  FOR  USE  TO  PRINT  CAPTION  IN  PLOT 

IHJMrTT 

EX=EXPt-TTe?HI) 

C  TEMPERATURE  =  TEMPS  IZP1)  AT  DEPTH  2=DEPIH  Z7 

TEM  =  MPS(ZPl) 

CALL  PLANCK?  TEM*  XV*  2*  PLNK) 

IF (.NOT.  DRUG)  GO  TO  31 

«RTVE<0.26!TT*X,AA.PHI,ex,XV.XLAH,DviN,PLNKtPLNKl,XWAVEN0 
WRITE (Q*32) TF«PS(ZP1 > ,Z7 

32  FORKAT?/"  TEMPS (7P1 > =«F10 . 3/»  Z7  DEPTH  AT  Z=*F10.3) 

31  CONTINUE 

C  NORMALIZE *  SEE  COMMENTS  IN  TERM!  ABOVE 

PLNK=PLNK/PLNK1 
T£RM2=PLNK»EX 

IF (.NOT.  OBUG)  GO  TO  33 
WRITE (0*341  PLNK.TERH2 
34  FORMAT (/•  PLNK  NORMALIZED  =»E12.S/ 

e  40X,  *  VALUE  OF  TE°M2  =«E12.S) 

33  CONTINUE 


NEXT  FINO  TEPM3 

TERM3  =  2«SUH  OF  FUNCTIONED  FOR  1=1 .2* . , . . ,M-1 
WHEN  DEPTH  INDEX  GOES  FROM  ZERO  TO  ?M*  MS  GOES  1  TO  M2P1 
THUS  TESH3  =  2«SUH  OF  FUNCTION (2J*1 )  FOR  3=1 » 2* . . . . .M-l 
SUMcO.O 
EXX=0.0 
MM1=H-1 
00  19  J=1.MM1 
KK=J 

M8=(2*J>*1 
I2TEMP=  JFIX (Z4) 

CALL  0EpTMP(«R*TEHPS*H2.I7TEHP.T£HPP,?CM} 

ZCH=ZCH/ZA 

CALL  THICK (TT*F «OS*ZCM,OENS> 

EX=EXP<-TT'PHI1 

C  SAVE  THE  FIRST  EX  IN  EXX  FOR  COMPARISON  TO  LATER  VALUES 

CALL  PLANCK (TEHPP*  XV,  2,  PLNK ) 

IF (.NOT,  OBUG)  GO  TO  35 

WRITE (0,26) TT*X,AA,PHI.EX*XV,XLAM*0WN»PLNK»PLNK1*XHSVEN0 
WRITE (0«36>  TcMPP* ZCh 

36  FORMAT  (/*  TEMPP=eF10.3/*  ZCM  DEPTH  =*F10.3 
3S  CONTINUE 

PLNKsPLNK/PLNK! 

$UM=SUM. (PLNK»EX) 


si? 

514 

515 
514 
S  1  7 

514 

519 

520 

521 
52? 
523 
S?4 

525 

526 

527 
325 
329 
53(i 
531 
53? 

533 

534 

535 
534 
537 
530 
539 
54  f. 
54* 

542 

543 

544 

545 

546 

547 
54? 

549 

550 

551 
55s 
55 » 
554 
55" 

556 

557 
550 
55° 
56  fl 
561 
56? 

563 

564 
=  65 

566 

567 


IF (.NOT.  OBUG)  GO  TO  3? 
WR!TE(Q,38)  PLNK.SUM 


n  non 
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38 

FORMAT  t/*  PLNK  NORMALIZED  =*E12.5/»  VALUE  OF  SUM 

INSIDE 

DEPTH  LO 

56B 

«0P  =*E 12.5  ) 

569 

3? 

CONTINUE 

570 

19 

CONTINUE 

571 

57? 

64 

CONTINUE 

573 

IF (.NOT.  DBU6)  30  TO  67 

5T4 

V.9JTefQ,6&>  KK 

575 

66 

FORMAT (*  THE  NUMBER  OF  TIMES  DEPTH  LOOP  COMPLETED 

BEFORE 

EXIT  IS  *14) 

576 

**K-> 

577 

6T 

CONTINUE 

578 

T£RH3=2.*SUH 

579 

560 

IF (.NOT.  DBUG)  60  TO  39 

581 

VRITE(Qo40»  TERM3 

58? 

40 

FORMAT (/ 

583 

<  SOX,  *  TERW3  =*£12.5  ) 

584 

39 

CONTINUE 

-585 

586 

587 

NEXT  FIND  TERM4 

588 

TERH4  =  4*SUK  Of  FUNCTION (21-1)  FOR  1*1 ,2, ,  .....M 

589 

SHEW  DEPTH  INDEX  COES  FROM  ZERO  TO  2H»  M.8  GOES  1 

TO  M2P1 

590 

THUS  TERH4  =  4*SL‘H  OF  FUNCTI0NC2J1  FOR  J=1 ,2, . e. . ,K 

591 

59? 

SUH=0.0 

593 

CXX=0.  C 

594 

00  20  J*t  tH 

595 

XK~J 

596 

M8=2*U 

597 

IZTEKP=  IFIXJ24) 

598 

CSti.  DEPTNP (K8 , TEMPS -M2, 1 ITEMP , TEMPR , ZCM ) 

599 

ZCM=ZCM/ZA 

600 

CALL  THJCK(TT,F,OS,ZCM,DtNS) 

601 

EX=EXPl-TT°PHn 

60? 

SAVE  THE  FIRST  EX  IN  EXX  FOR  COMPARISON  TO  LATER 

VALUES 

603 

CALL  PLANCK (TEKPP, XV,  2,  PLNK) 

604 

60H 

IF  t .NOT .  D8UG)  GO  TO  41 

606 

S»RITE{Q,26)TTtX,AA,PHt,EX,XV?XLAM,0WN,PLNK,PLNKl*XMAVEN0 

607 

VRJTE(0,36)  TEMPR, ZCH 

608 

41 

CONTINUE 

609 

610 

PLNK=PLNK/PLNK1 

611 

SUKsSUH, <PLNK«EX) 

612 

613 

IF(oN9T.  DBUGJ  GO  TO  42 

614 

WRITE (0,38) PLNK, SUM 

615 

42 

CONTINUE 

616 

617 

20 

CONTINUE 

618 

65 

CONTINUE 

619 

IF (.NOT.  05UG)  GO  TO  68 

620 

VRITE(G^66>  KK 

621 

88 

CONTINUE 

622 

non  non  no  o  no  non 


TERH4=4<,*SUM  623 

624 

IF  (.NOT.  D8UG1  GO  TO  43  625 

URITE!0.44)  TERM4  626 

44  FORMAT {/  62? 

*  60X*  '»  TERM4  ="£12,5  J  628 

43  CONTINUE  629 

630 

631 

NEXT  FIND  INTEGRAL  632 

INTEGRAL  =  <H/3)*(TERM1  *  TERN2  ♦  TERH3  •»  TERM4  )  633 

H=(B-4I/2M  63* 

H=(Z-0«)/FLOAT(P2>  635 

IN  CUR  PROBLEM*  ME  ME  SET  BOUNDARY  A  TO  ZERO.  636 

SUM  =  INTEGRAL  63? 

5UM=0.0  638 

SUK=<H  /3. 1*  <T£RM1  ♦  TERK2  ♦  TERH3  ♦  TERM4  )  639 

EMERGENT  INTENSITY  EINT  FOR  FREQUENCY  IN  THIS  LOOP  FOLLOWS  640 

EINT <JJJ}  a  PHI*SUH  641 

642 

IFt.NOT.  DBUG)  GO  TO  45  643 

WRITE (0.46) TERM! • T£RH2 »  TERH3 *  T£RH4 , SUH  t  PHI *  E INT ( J J J )  644 

'.6  FORMATS/*  TERH1  =«El?.5/»  TERM2  s*EI2,5/»  TERM3  =*Ei2.S/  645 

*  *  TERH4  =*E12.5/*  SUH  THE  INTEGRAL  »«E12.S/  646 

*  *  PHI  »*EI2.5  /*  EMERGENT  INTENSITY  =*E12.5*/////  )  647 

45  CONTINUE  64 6 

649 

>99  CONTINUE  650 

5SE  HAVE  COMPLETED  COMPUTATION  FOR  A  GIVEN  FREQUENCY  651 

NEXT  WE  INDEX  FREQUENCY  DO.  LOOP  AND  REPEAT  652 

653 

654 

WRITE  THE  OUTPUT  655 

KY3  CONTAINS  WAVELENGTH  IN  ANGSTROMS.  656 

XVV  IS  EQUIVALENCED  WITH  XV3.  657 

WRITEtO.U)  658 

659 

PLNX2  IS  PLANCK  VALUE  AT  CENTER  FREQUENCY  XL AH  AND  TEMPERATURE  660 

TEMPS (ZHP11  TO  WHICH  MAXIMUM  OF  EMERGENT  POWER  IS  NORMALIZED.  661 

LATER.  SEE  BELOW.  662 

TEK=T£MPS (ZMP1 #  663 

CALL  PLANCK {  TEK.XLAM,  2.  PLNK2)  664 

665 

C  INTEGRATE  THE  SPECTRUM  AND  WRITE  THE  YALUE.  666 

VAL=0.  667 

NFQHI=NFREO-l  668 

DO  95  J=1.NFGM1  669 

95  VAL=(((EINT(J)<-FlNT{J»l))/2.)*(XVV(J*l)-XVyUn>*VAL  670 

WRITE(Q*96)  XVV I i ) *XVV (NFREG) »  VAL  671 

96  FORMAT!/*  THE  EMERGENT  POWER  INTEGRAl  ./UK  NON-NORHALIZED)  672 

•AS  ORIGINALLY  COMPUTED*  FROM*  673 

*  /H2.3*  ANGSTROMS  TO*  F12.3*  ANGSTROMS  IS*  EI2.5)  674 

WRITE(Q*102)  675 

102  FORMAT <*  POWEP  VALUES  AS  COMPUTED.  NOT  NORMAL IZ£0»)  676 

WRITE  t-3*  85  i  677 


on  no  no  on 
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85  FORHff!//*  EMERGENT  INTENSITY*)  678 

WRITE  <Q* 47) (!<XyV!J) »EINT<J) ) ,J=1 .NFREO) )  679 

67  FORMAT !//{5!2x*F12.3,F12.5)?)  680 

681 

PLOT  OUTPUT  68? 

NORMALIZE  MAXIMUM  OF  EINT  TO  UNITY  SO  WE  CAN  PLOT  683 

YMAX=EINT { i )  684 

DO  56  J=2. NFREO  685 

IF  (EINT  { J)  .GT .  YMAX)  YHAX=EINT ! J)  686 

58  CONTINUE  637 

NOW  YMAX  CONTAINS  THE  NORMALIZATION  CONSTANT  FOR  EINT  688 

689 

NEXT  FIN05  NCAROS  The  NUMBER  OF  DATA  CARDS  IN  A  SET.  690 

AKKH=FLOAT INFREO)  691 

AKKKH=AKXH/4.  692 

NFRE04=NFRE0/4  693 

I F ( ASS ! AKKKH-NFRE04 J  .LT.  0.2)  104,105  694 

104  NCARDS=NFREQ4  695 

GO  TO  106  696 

105  NCARDS=NFRE04 ♦ 1  697 

106  CONTINUE  698 

699 

Al  IS  AA  TIMES  T0RR/760  700 

A1=!AA*TORR)/?60.  701 

702 

NEXT  DO  LOOP  PUNCHES  NOUPS  SETS  OF  CARDS  703 

IF1PUNCM1 ) 107 » 1 08  704 

1C7  CONTINUE  705 

00  103  KA=1 ,NDUPS  706 

JT1 =4  707 

J8B=1  708 

PUNCH  109, DM, YEtZtDENS.Al, GROUP, TORR  709 

109  FORMAT (lX,£'10, A2«*Z=*F3. 1 ,*DENS=*E12.5,*  AA=*F7.5,*  GROUP*I3,  710 

*  F6.0,»T0RR  RADIANT*)  711 

PUNCH  110, VAL*  XVV (1 ) ,XVV INFREO) «  YK4XyTEH,PLNK2,NFREQ,NCARDS  712 

110  FORMAT (lx,E12.5«2F9»3,E12.5,F6.0*E12»5,2I5  >  713 

INTEGRAL  VALUE  HERE  IS  NOT  NORMALIZED.  714 

DATA  HERE  ARE  AS  COMPUTED  715 

00  111  L=l» NCAROS  716 

PUNCH  112, L,  !<XVV!J),EINT(J)>,J=JBB.JTT>  717 

112  FORMAT!  14 .4 (F7,2*E3 2.5) )  718 

J88=JB8‘4  719 

JTT=JTT  *4  720 

111  CONTINUE  721 

PUNCH  113,  22  722 

113  FORMAT  !A1) 

103  CONTINUE  724 

108  CONTINUE  725 

726 

C  NEXT  00  THE  NORMALIZATION  727 

DO  61  J=I»NFRF.O  728 

EINT(J)=EINT!J)/YMAX  729 

61  C0NTINUE  730 

IF!  .NOT.  PLOT!  )  GO  TO  57  731 

C  EINT  IS  NORMALIZED  ANO  IS  READY  FOR  PLOT  73? 


no  a  r>  non 
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WRITE(O.ll)  733 

WRITE (0.59)  734 

59  FORMAT UH1.50X,  «  FROM  PROGRAM  LTE4  *)  735 

WRITE (Q. 63)  DM.YE. IHR.MIN. I SEC  736 

*  .GROUP'.  TORR  737 

WRITE (0.24) H.Z.F.0S.DCN5* AA.NFREQ.TEMPL.DVN.XLAH.D8U6.STEP.PL0T]  735 

WRIT£(O.B1)OS1 .0S2.XLAM] .XLAM2  739 

WRIT£<0.89)  TEMPS ( 1 ) .  TEMPS!2HP15  740 

89  FORMAT (*  BOUNDARY  TEMPERATURE  tS»F7.0»3X,e  CENTER  TEMPERATURE  IS*  74] 

*  F7.0)  74? 

WRITE(O.IO)  062.(TEKPS(J),J=1,ZP1)  743 

WRITE(Q.22)  ft  DWN.AHW  744 

WRIT£(Q.23)  XLAM.PLNK1  ’45 

PLN3C2  IS  PLANCK  VALUE  AT  CENTER  FREOUENCY  XL AM  ANO  TEMPERATURE  746 

TEMPS (ZHP1J  TO  WHICH  MAXIMUM  OF  EMERGENT  POWER  IS  NORMALIZED.  747 

SEE  BELOW.  748 

WRIT£(G.98)  PLNX2.  XLAM.  TEM  749 

98  FCRHAT(«  PLANCK  VALUE  IS*  £12.5*  AT  CENTER  FREOUENCY*  F12.3/  750 

*  *  TND  AT  CENTER  TEMPERATURE*  F12.3  'DEGREES  KELVIN.*/  751 

*»  TO  WHICH  MAXIMUM  OF  ENE-.ftENT  POWER  IS  NORMALIZED.*  )  75? 

WRITE(0.60)  YHAX.DUH  753 

6C  FORMAT («  YHAX  THE  NORMALIZATION  CONSTANT  FOR  EINT  IS*EI2.5./  754 

*  *  THE  OPTICAL  THICKNESS  (TOTAL)  IS*E12.5>  755 

CALL  PLOT (EINT.XVV.NFREQ)  756 

57  CONTINUE  757 

758 

NEXT  NORMALIZE  TO  PLNK2  VALUE  759 

DO  97  J=1.NFRE0  760 

INTEGRATE  NEXT  76? 

97  EINT ( J) =  EINT ( J)  «  PLNK2  76] 

VAL=0.  763 

NFQHI=NFRE0-1  764 

DO  100  J=I .NFOM1  765 

.90  VAL=(<(EINT(J)*EINTOI)  ) /2. ) *<XV¥ ( J+I J-XVV < J> > > *VAL  766 

WRITE(Q.Il)  767 

WRITE(Q»98)  PLNK?.  XLfiH,  TEH  768 

WRITE(Q.99>  XVV( I ) *XVy (NFREQ) .  VAL  769 

99  FORMAT (  /*  TH *  EMERGENT  POWER  (MAXIMUH  NORMALIZED  TO  PLANCK)  FRO  770 

*0M*  F12.3* ANGSTf  OMS  TO  *FJ2.3*  ANGSTROMS*/  771 

*  «  IS*EI 2.5  )  77? 

WRITE (0. 101 )  773 

101  FORMAT (»  POWER  NORMALIZED  TO  PLANCK  VALUE  FOLLOWS*)  774 

WRITE(0.8S)  775 

WRITE(G. 47) (((XVV(J),EINT(J)),J=1. NFREQ))  776 

777 

PUNCH2  NEXT  ,'78 

NEXT  DO  l OOP  PUNCHES  NDUPS  SETS  OF  CARDS  779 

IF (PUSCH2) 114.115  780 

114  CONTINUE  78] 

DO  116  KA-1. NDUPS  73? 

J7T=4  703 

J5P=1  784 

PUNCH  I09.DP.YE.Z.0ENS.A1 .GROUP. TORR  785 

PUNC't  XVV(l)  .XVV(NFREO)  »  YM AX ♦  TEH. PLNK2. NFREQ »NC ARDS  786 

C  INTEGRAL  VALUE  IS  AFTER  NORMALIZATION  TO  PLNK2.  787 
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C  DATA  ARE  NORMALIZED  TO  PLNK2.  788 

DO  117  L=1.NCARDS  789 

PUNCH  112, L,  ( ( XVV ( J)  *EINT ( J) >  *  >I=JBB»  JTT)  79P 

J8B=J8B*4  791 

JTT=JTT*4  792 

117  CONTINUE  793 

PUNCH  113,  ZZ  794 

116  CONTINUE  795 

115  CONTINUE  796 

797 

C  REUSE  XVI  STORAGE  TO  HOLD  WAVELENGTHS  FROM  EYEBAL  798 

CALL  EYEOAL  (XVI, 392)  799 

800 

C  NOW  PUT  A2  IN  C0MM0N/B6/  INTO  PHI1  STORAGE  801 

00  119  J=1 *392  '  802 

PHI  1  ( J)  =  A2  (-J5  803 

119  CONTINUE  804 

BOS 

CALL  COREKT (EINT*XVV,NF REQfPHIl ,XV1 ,392,1 , PHI 2)  806 

C  B07 

C  PUT  CORRECTED  DATA  FROM  PHI2  INTO  EINT  808 

DO  120  J=1*NFR£Q  809 

EINT<J1=PHI2(JJ  310 

120  CONTINUE  811 

C  812 

C  INTEGRATE  THE  LUMINOUS  VALUE  AND  WRITE  THE  VALUE  813 

VAL=0.  '014 

NF0M1=NFRE0-1  815 

DO  129  J=l,NFOMl  816 

129  VAL=!HEIN7(J).EINT(J+in/2.)»(XVV(J*l>-XVV<J))»+VAL  817 

WRITE (0, 130 > XVVO  5  *XVV(NFREQ) »  VAL  813 

130  FORMAT </»  THE  LUMINOUS  POWER  INTEGRA . (MAXIMUM  NON-NORNALIZED)  819 

•AS  ORIGINALLY  COMPUTED*  FROM6  820 

•  ,'*712.3«  ANGSTROMS  T0«  F12.3*  ANGSTROMS  ’S*  E12.5)  821 

822 

C  FIND  H»X  OF  EINT  LUMINOUS  ENERGY  823 

YMAX=EINT (1 J  824 

00  12T  J=2*NFREQ  825 

IF (EINT ( J)  *GT.  YMAX)  YMAX=EINT ( J)  826 

121  CONTINUE  827 

C  NOW  YMAX  CONTAINS  NORMALIZATION  CONSTANT  OF  LUMINOUS  POWER  EINT  828 

C  8  *9 

WRITE (0* 123)  830 

123  FORMAT!///*  RELATIVE  LUMINOUS  POWER  NOT  NORMALIZED*)  831 

WRITE(Q,47) ( < (XVV( J) »EINT ( J) ) ♦ J=1 tNFREQ) )  832 

C  833 

C  NEXT  NORMALIZE  EINT  TO  UNITY  834 

DO  122  JslcNFREO  835 

EINT (J)=EINT ( J) /YMAX  836 

122  CONTINUE  837 

C  838 

C  PUNCH3  HERE  839 

C  TH'  S  IS  LUMINOUS  POW.L  •  JF.MALIZED  TO  UNITY  840 

C  NEXT  00  LOL  P  PUNCHES  NUUPS  SETS  OF  CARDS  841 

IF (PUNCH3) 125,126  84? 


125  CONTINUE  043 

DO  127  KA=1.NDUPS  044 

JTT=4  045 

J8B=1  846 

PUNCH  124.0M.YE.Z.0ENS*  Al  .GRO'JP.TORR  847 

124  rORMATUX.AiO.A2.*Z=*F3.1»»DENS=«E12.S.»  AA=>*F7.5.*  GR0UP»I3,  048 

*  F6.0.*7QRR  LUMINOUS*)  849 

PUNCH  110.VAL*  XVVtl>,XVV{NFREO>t  TMAX.TEH.PLN<2*NFREQ»NCAROS  850 

C  INTEGRAL  VALUE  VAL  IS  LUMINOUS  BEFORE  NORMALIZATION.  851 

00  128  *.=1  »NCAROS  852 

PUNCH  112. L«  ( (XVV<  .EINT ( J) )  *  JaJ8S»  JT''5  853 

JB8=JBB*4  854 

JTT=JTT*4  855 

128  CONTINUE  856 

PUNCH  113.  ZZ  857 

127  CONTINUE  858 

126  CONTINUE  859 

WRITE (Q.49)  860 

49  FORMAT tlHl .10 </) .SOX**  END  OF  CASE*  1  861 

GO  TO  1  862 

END  863 


n  n  oo  oo  on 
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SU8R0UTINE  PLOT (EINT »XVV»NFREQ)  864 

C  WRITTEN  BY  B  E  DOUDA  4  SEPT  1969  86* 

868 

DIMENSION  LIN£<132)»EINT(2000)  «XVV (2000 >  869 

INTEGER  0  870 

DATA  Q/6/  8  T1 

DATA(8LANK=IH  >»  (D0L=1HS),  (X=1HX),  {  Y=1H») *  (  Z=1H0>  87? 

TYPE  INTEGER  BLANK,  00L,  X,Y«  Z  8f3 

C  874 

WRIT£{0«9)  875 

9  FORMAT </45X38HRELATIVE  POWER  (NORMALIZED)  876 

*/,16X3H0.022X4H0.252IX3H0.522X4H0.752SX3H1.0  )  877 

878 

PRINT  A  LINE  OF  OOL  TO  H4KE  VERTICAL  AXIS  879 

00  10  J=1 , 3.32  880 

10  LINE ( J) =8LANK  B81 

00  20  J=16, 1 19  88? 

20  LINE(J>  =  OOL  883 

WRITE<0.30)  (LINE(J),J=16,H9)  884 

30  FORMAT  (4X9HANGSTR0MS2X104A1)  885 

886 

BLANK  THE  LINE  887 

00  40  J=17.1ie  888 

40  LINE(J)  =  BLANK  889 

NFRsNFREQ  890 

DO  70  N=1.NFR  891 

892 

PUT  X  IN  SELECTED  LOCATION  893 

IF <  EINT (N)  .LT.  0.  .OR.  EINT (N)  .GT.  1.  )  GO  TO  120  894 

J=I00.  *  EINT (N)  ♦  18. S  R9S 

LINE(J)  -  X  896 

SO  CONTINUE  097 

WRITE(0»60 )  N»XVV(NJ»(LINEIL)*L=i6,119)  898 

60  FORMAT (I5,1X,F8.2»1X* 104A1 >  e99 

900 

PUT  A  BLANK  BACK  IN  THE  LOCATION  WHERE  X  WAS  901 

LINE(J)  =  8LANK  902 

70  CONTINUE  903 

DO  SO  J=l«132  904 

80  LINE(J)=BLANK  90S 

DO  90  J=16* 1 19  906 

90  LINE(J)  =  DOL  907 

C  908 

WRITE(Q.IOO) (  LINE(J).J316.119)  909 

100  FORMAT (15X10461)  910 

Pr'URN  911 

c  91? 

120  fe^ITEIO.nO)  913 

130  FORMAT (//IX, *EINT  IS  NEGATIVE  OR  GREATER  THAN  1.0.  IF  NEGATIVE  914 
•.  LOOK  FOR  ERROR  IN  PROGRAM.*?*  IF  GREATER  THAN  ONE,  915 

*  CHECK  NORMALIZATION*//)  916 

RETURN  917 

END  918 
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SU8R0UTINE  PLANCK  <T»X»N,R)  919 

C  PROGRAMMED  BY  B  E  OOUDA.  OCT  1970.  RUNS  ON  CDC6600.  920 

C  T=TEHPERATUR£,  X=HAVELENGTH  IN  ANGSTROMS,  R=RADIANT  EXITANCE  921 

C  N  DETERMINES  UNITS  ASSIGNED  TO  R  92? 

C  C1=3.7415E-S  CH*»2  ERG/SEC  =  3.741SE11  A**2  ERG/SEC  923 

C  C1-3.7415E4  A»*2  J/SEC  =  H  A«»2  ci=2»PI*H*C«C»  924 

C  C2=l. 43879  CM  K  =  1.43879E8  A  K  925 

C  R  =  WATTS/A**?/A  R  VALUES  ARE  PER  HEMISPHERE.  Q26 

C  R*1.E20  =  WATTS/CM**2/K,CR0N  WHEN  N=2  92T 

C  THERE  ARE  IE*8  ANGSTROMS/CH  928 

C  ONE  MICRON  =  IE+4  ANGSTROMS  929 

C  R  =  PLANCK  POWER  930 

DATA  (Cl=3.741SE4>c  (C2=I .43P79E81  931 

E  =  EXP {  C2/fX*TU  ?>3r 

P-=  Cl'f  <.V.**S.)  *  '£-!,) X  933 

IF (N  .EO.  2  )'?=R»1.E20  93'. 

RETURN  935 

END  93fi 
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SUBROUTINE  CEPTMP1H8. TEMPS tR2»Z.TE«PP.ZCM) 

937 

938 

c 

B  E  DOUDA  25MA°CH  1972 

939 

940 

OIMt  ' 3JON  TEMPS (201 ) 

941 

INTEGER  7 

94? 

INTEGER  7P1 

943 

944 

c 

FIND  DEPTH  CALLED  SIZE  FCR 

EACH  SIMPSON 

INTERVAL 

94S 

SIZE=FL0AT(Z)/FL0AT(M2) 

946 

ZCH=(M8-I!  *  SIZE 

947 

c 

ZCH*  THE  DEPTH  IN  CM  AT  H8* 

IS  RETURNED 

TO  MAIN  PROGRAM 

948 

949 

c 

NEXT  FIND  TEMPP  AT  H8 

950 

ZPI=Z*1 

951 

M2P1=M2*1 

95? 

c 

HERE  SIZE  IS  THE  NUMBER  OF 

SIMPSON  RULE 

INTERVALS/CM  DEPTH 

953 

SIZ£=FLOAT (N2) /FLOAT !Z) 

954 

IF(H8  *EQ.  1 )  GO  TO  30 

955 

IF(H8  ,EQ.  H2P1)  GO  TO  40 

956 

00  10  J=2.ZP1 

957 

K-J 

958 

X=(J-1»*SIZE 

959 

IF<X  „GT.  (H8-11)  GO  TO  20 

960 

iO 

CONTINUE 

961 

20 

DIFF=X»(M8-1) 

96? 

FRAC=D IFF/SIZE 

9S3 

DIFF=TEHPS(K)  -  TEMPS(K-l) 

964 

SIZE=FRAC*DIFF 

965 

TENPP  =  TEMPS(K)  -  SIZE 

966 

RETURN 

967 

968 

30 

TEMPP  a  TEMPS ( I) 

969 

RETURN 

970 

971 

40 

TEMPP  =  TEMPS (ZPI  ) 

972 

RETURN 

973 

END 

974 

finooo 
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function  VOIGT (X»A) 

PROGRAM  WRITTEN  BY  G.  RY8ICK I 

COMPUTES  VOIGT  FUNCTION  FOR  ANY  VALUE  OF  X  AND  ANY 
POSITIVE  VALUE  OF  A.  BEFORE  CALLING  FIRST  I 3HE 
ENTER  COFVOI  WITH  THE  STATEMENT 
DUMMY! =COFVOI <DUHMY2*DUHKY3) 

DIMENSION  C C31 > 

TYPE  COMPLEX  Z 
IF(A.EQ.O.O)GO  TO  50 
Al=3.«A 
A2=A*A 

IFJA.LT.O.DGO  TO  10 

Z=CEXP (-9,424777960769360  A* 9. 424777960?6938*X) 

VOIGT=0,0 

GO  TO  20 

10  Z=CCOS (9.42477796076938*X*9.42477796C7693B*A) 
VOIGT=0.564189583547756*CXP(A2-X*X)*COS(2.*A*XJ 
20  81= I I .-REAL (2))*A*1.5 
B2=-A!KAG(Z) 

S=-0.-1.5*X 

T=S«S»2.2S«A2 

DO  40  N=1 »31 

T=T*S*0.25 

S=S»0.5 

B1*A1-BI 

B2=-B2 

IF (T*GT .2.5E-I2) GO  TO  30 
VOIGT=VOIGT-C(N)*A«. 333333333333333 
GO  TO  40 

30  VGIGT=VOIGT*C (N) *(81 ♦B2*SJ /T 
*0  CONTINUE 
RETURN 

50  VOIGT=0.S64I89583547756«EXP (-X*X> 
return 

ENTRY  COFVOI 
K=-16 

00  60  H=I*31 
K=K*1 

60  C(NJ =0„0897935610625833»EXP (-FLOAT (K*KJ/9,j 
RETURN 
END 


975 

976 

977 
976 

979 

980 

981 

982 

983 

984 

985 

986 
98? 

988 

989 

990 

991 

992 

993 

994 

995 

996 
99? 
993 
999 

1000 

10C1 

1002 

1003 

1004 

1005 

1006 
1007 
1006 

1009 

1010 
ion 
1012 

1013 

1014 

1015 


c 

c 

c 


c 


SUBROUTINE  THICK <TT*c,05*Z*Dt:NS> 

INPUT  OS*  Z-  OEMS 

OUTPUT  TOTAL  OPTICAL  TMlC.NESS  . f  AT  DEPTH  2 
B.  E.  DOUDA  20  MARCH  72 

DATA (P=3. 141592654) .  <E--4.30298E-010)  *  (HM=9. 1091E-C2B) 
OATA(C’'2.99?925E*0I<>> 

DATA <7TC=6.00CE* 007) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

C 

c 

C 

c 

c 


TT=  TOTAL  OPTICAL  ThICKNESS=  DIMENSIONLESS 

P=  PI 


E=  ELECTROSTATIC  CHASGE  {CH*#1.5>  (S»*-l> 

Z-  PHYSICAL  CEPTH  Or  ATMOSPHERE  CM 
P=  DOPPLER  HALE  WIDTH  IN  FREQUENCE  UN'ITS 
RM=  ELECTRON  REST  MASS  G 

SPEED  OF  LIGHT  <CK>  <S**-1) 

DEMS=  N12-*  PAPTtCLES/CC  IN  ATOKSPHERE 

OS=  OSCILLATOR  STRENGTH  CF  LINE*  DIMENSIONLESS 

H=  PLANCK  CONSTANT  ERG  S 

V=  CENTER  FREQUENCY  OF  LINE  1/S 

ERG=  DYNE  CM  DYNE=  GH  C?’'SEC**2 


B12»  EINSTEIN  COEFFICIENT  =■  <OS«4*P«P«l£«£)/<F«RM»C) 
ABOVE  B12  IS  FORM  WHEN  N12  IS  IK  DENSITY  UNITS 

TT=  <H«V)  /  <4*P*F>  *  (N12  *  fil2  -  N2l  *  021)  «  , 
<21  *  B2I  IS  NEGLIGIBLE 


TI  a  <OENS*OS  *P  *  E  *  £  *  Z  >  /  <  F  *  RM  «  C  ) 


C  OEySC=  DENSITY  AT  SPECIFIED  THICKNESS  ,  TTC 

C  DEI.  ;  a  (TTC*F*RP*C)/<OS  *P«E«E«Z) 

C  WRITE  <61,101  OENSC*  TTC 

C  10  FOPHAT  fiHi,2E21.14) 

RETU2N 

FND 


1016 

1017 

lOlB 

1019 

1020 
1021 
1022 

1023 

1024 

1025 

1026 

1027 

1028 

1029 

1030 

1031 

1032 

1033 

1034 

1035 

1036 

1037 
103R 

1039 

1040 

1041 

1042 

1043 

1044 

1045 

1046 

1047 

1048 

1049 

1050 

1051 

1052 
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SUBR0UT1NE  OOPPLERIXLAH.NK.TEKP.PHWIDTHtDWN  ,AHW  ) 

C  COMPUTES  DOPPLER  HALFWIDTh  OF  ATOMIC  LINE  AT  SOME  TEMPERATURE 
C  WRITTEN  BY  f  E  QOUOA  25  AUG  1969 

C  TEMP  IS  DEGREES  KELVIN*  XHASS  IS  WEIGHT  PER  MOLE* AND 

C  XlAMBOA  is  THE  CENTER  OP  THE  RESONANCE  LINE  IN  ANGSTROMS. 

r 

V 

OIHENSION  AHG<  1 « 1 ) »XLAHBDA(2) < AHHIOTHO) <0WAVENO (3» 

DIMENSION  FP.EO (  1) 

TYPE  REAL  LAHBDAPt  LAMBDAM 
'S'  INK  .WE.  0  )  GO  TO  10 
XKASS  2.9898 

C  XLAHBDA<2>  a  5895.9236 

C  XLAM3DA ( 1 )  =  5389,9504 

L'  XLAMBDAII  1=5893. 

C  XLAMBDA(2)=S893. 

K=l 

XLAUBDA (K) =XLAH 

C  XLAN  =  LINE  CENTER  IN  ANGSTROMS, 

a0LT2-1.38054E-0l6 
SPEECH  =2.9979?5EO10 
AVAGAD=6 . 022S2EC23 

CONSTAN  =2.0*S0RT  <2.0*80LTZ)*SQRT ( ALOG (2.0 1 ) ‘SORT ( AVAGAO) /SPEEDLI 

C  FHUIDTH  15  DOPPLER  HALFWIDTH  IN  FREQUENCY  UNITS 

FKWIDTH= (CONSTaN-SPEEOLI/ (XLAMBDA (K>  *1 .OE-8) J  *SQRT (TEHP/XMASS) 
XWAVENO  =1.0/(XLAHBOA(K)  •  I. OE-8 J 

DWAVENO(K) =FHWIDTH/SPEEDLI 
C  OWN  a  DOPPLER  WIDTH  IN  WAVENUMBERS. 

OWN=OWAVENO(K) 

HWAV£NO=  0WAVEN0(K>/2.3 
LAHSDA?=1.0E«8/  (XWAVENO  ♦  HWAVENO) 

LAH&DAH=l.GE*8/r  (XWAVENO  <»  HWAVENO) 

C 

C  AHWICfH  IS  DOPPLER  HALFWIDTH  IN  ANGSTROMS 

A'  WIDTH(K)  a  LAMBDAM  -  LAH8DAP 
Mb'  =  AHWIDTH(K) 

RETURN 
10  CONTINUE 
NFI=HFREO-I 
DO  cSS  1=1  »NFR£G 
TO=NFftta-I*I 

ANG(I*NF1*K)  =  1.0E*8  /  (XWAVENO  -  (OWAVENO(K)  *  FREQ(I))) 
ANGUD.K)  =  1.0E«8  /  (XWAVENO  ♦  (DWAVENOOO  t  FREO(I)!) 

20  CONTINUE 

AHWIDTH (3) k^HWIDTH (2> 

J=K 


1053 

1054 

1055 

1056 
1059 
J-060 
1061 
1062 

1063 

1064 

1065 

1066 

1067 

1068 

1069 

1070 

1071 

1072 

1073 

1074 

1075 

1076 

1077 

1078 

1079 

1080 
1081 
1082 

1083 

1084 

1085 

1086 

1087 

1088 

1089 

1090 

1091 

1092 

1093 

1094 

1095 
TG96 

1097 

1098 

1099 

1100 
1101 
1102 

1103 

1104 

1105 

1106 

1107 

1108 


OWAVEtfOO)  =  DWAVcNO(2) 
C  CALL  PRINT (4) 

RETURN 

END 


I10« 

1110 

nn 

1112 

111? 


onoooooooo  non 
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subroutine  COREXT  (A.XA.KA.B.XB.KBtH.C)  1U4 

PROGRAMMED  by  a  c  DOUDA.  OCT  1P70.  RUNS  ON  CDC6600.  Ills 

DIMENSION  AUOOOHXAUOOQ)  ♦  B ( 1000) . XB( 1000) ,  CC1000)  1116 

DIMENSION  A(4000) .XA<4000) .  8(4000) «X8(4000>  »  C<4000)  HIT 

DIMENSION  A (20 00) »XA (2000) t  B<2000> «X9(200Q) ,  CI2000J  1118 

INPUT.  A(XA)  ARE  INTENSITIES  TO  BE  CORRECTED  AT  KA  POINTS  1U*> 

INPUT.  XA(KA)  ARE  WAVELENGTHS  ASSIGNED  TO  A  H20 

INPUT.  B(K8)  ARE  DATA  THAT  CORRECT  A.  1321 

INPUT.  XB (KB)  ARE  WAVELENGTHS  ASSIGNED  TO  B  1122 

INPUT.  KA  IS  NUHBER  OF  A,  X8  IS  THE  NUMBER  OF  B  1123 

H  SETS  TYPE  OF  CORRECTION  112* 

DOES  A/B=C  AT  WAVELENGTHS  ASSIGNED  TO  A  WHEN  H=0  1125 

DOES  A«B=C  AT  WAVELENGTHS  ASSIGNED  TO  A  WHEN  K=1  112ft 

OUTPUT.  C(KA)  ARE  THE  CORRECTED  DATA  AT  KA  POINTS  AND  1127 

AT  WAVELENGTH  XA.  ll2S 

IFCXA(H3LT.XBai.GR.XA(KA>,GT.XB(KB>)  WRITE(6.1)  1*29 

1  FORMAT*//. *  CHECK  RANGE  OF  B  IN  SUBROUTINE  COREKT  FOR  ADEQUACY*//)  U30 
X8M=KB-1  1131 

DO  3  JH.KA  1132 

DO  3  K=1«KBM  H33 

KP=K*1  1134 

IF!  .NOT * (XB(K) .LE.XA ( J)  .A'ND.X3*XP)  .GF.XAIJ)))  GO  TO  3  H3S 

D=XB(KP)-Xfl(K)  H36 

E=XA(J>-X3(X)  H37 

F=B (K) ♦< (E/D) ®(B(KP)-8(K)))  1138 

IF (M.EO.l >  GO  TO  5  1139 

C(J)=6(J)/F  H40 

GO  TO  3  1X41 

5  C  ( J)  =A  ( J)  <*F  1142 

3  CONTINUE  1143 

RETURN  1144 

END  1145 
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8L OCK  DATA  EYES  U*6 

COMMON/36/ftl {392)  1)47 

C  Ai  t S  EYE  RESPONSE  TABLE  USED  IN  EY£8AU  114a 

OATAUAllJ), 3=1,1021  =  1149 

*0. 00004000  ,0.00004499.6.00004947 .0. 00005 :>88*0.80f,a586e.0 .00906425*  H50 

*0,000071 10, 0.90007966 .0,06005037, 0.00010367*0. 00912000. 0.00013763,  1151 

«O„0CeI549S*0, 000)7290* 0,30019260*0.00021506*0-00024100* 0.63027150,  1152 

«0, 00035750, 0.00035000, 0.00040000, 0.00044990*0. 00049470, 0.00053380,  1153 

*0,00058660*0.00 164250*0. 00071 100, 0.06079660, 0.00090370,0. 00103670,  1154 

*0,00120030, 0.00237740,0. 00155520, 0.00174189, 0.00194560*0. 00217500,  1155 

*0.00243540, 0.00274420, 0.00310080, 9, 00351660, 0.004C00C0, 0, 00454050,  1156 

*6.00515200.0, QOSilOTSO.O^OOSSmO^,  00726250, 0.00305600, 6. 0088B<:S>J*  1137 

*0. 00976000, 0, QH664S0, 0,01160000, 0,01256990,0. 01357920,0, 01462730.  1158 

*0,01571360*0. 0168375-3, 0*01 799840, 0*01919570. 0.02042880.0.02169710,  1)59 

*O,O23eO0OC*3. 02432270*0, 02565760,0.02701490, 0,02840480, 0.02963753,  1160 

°0. 03132320, 0.03237210*0.03449640*0,0362063!'  *0.03806000,0. 03987600.  1161 

*0.04180800,0.04380200, 0,0452-640 0*0. 04800000 »O,05<'2! 600 ,0.05251 800*  1162 

*0. 0549124?, n. 05740400, 0.66000000,0.06265397, 0.065351 96, 0.06810296,  >163 

*0.07093596 *0.07387495 ,0.07694359,0 ,0301 8700 *0.03356795, 0.087) 7096,  ) ,64 

*3.89099996, 0.05501696, 0.00917599,0. 10349898,  >'.10796797, 0.1 1262494,  1165 

*6.11^471 99*11.3  2252598,0. 12773395  >0.13327295,0. 13900000,0.14494245)  1166 

DATA)  <Al 131  ,  J=1 03*20*,)  -  1167 

*0.35071944, 0. 15570747*0, 16287994, 0.16931248, 0.17607949,0. 18325746,  116?. 

•I), 19052500,0  1491 4246, 0.20 799994,0 .2 173# 446,0.22701597,0.23707145,  1159 

•'0.24756795. 0.25836245,0. 1701 1198, 0.28227347,0. 29510397, 0.30866045,  1|70 

*0.32299995, 0-33823997, 0.35440999, 0.37137997,0.38901997. 0.40724959,  1371 

*0.42593998,0.4445608., 0.46422994. 0.48361999, 0.50299996*0.52285999,  1172 

*0.54354995, 0,56463996.0.^649999, 0.60324996*0, 62985998,0.65109998,  1173 

*9. 6?ic8999,0.65n9999,«.  *.3599998, 0.72770000, 0.74487996,0.76354000.  1174 

*0.77763999*0. 7923 8994* &.808 I 5^95,0. 52253599 r 0,83631 998 *0-8494795*,  1375 

*0. 86199999*4. '7385994, 0.83505995,0. 8550 1999. 0,90556997, 0.91*94000,  3 176 

®C. 92374958,5, 53293995, 0.939619<>, 0.5471399?, 0,95400000, 0,96034998,  1177 

*0. V6S1195S, O.V?:3399a, 0,97603955. 0.96024994,0. 98399997,0. 98733 9?S,  1178 

*0. 99023999, 0,3927895§,0,9S4!^9!?a*t-.9966&39S*0.9582l20C,  0,59536795.  117“ 

*1. ODOCIS26. 1.00029966. 1, 00006390, 0. 99943199*0. 98-40?^6, 0.59692595*  1130 

*0.99499595  *  0 . 0926 1 599 , 0 , 55977995  *0,93648995*0, 98276097 , 0.9786299?*  11 S 1 

*0.97436995, 0.56912998, 9. 96377999, 0,9580799^*9, 95195996, 0.04S5J998,  3 182 

*0,93358999,0. 93124998,0.32369594. 0.91537994*0. 93589999, $ „89gl 0947,  l'$3 

*0.58900995 ,0.87963°98* 0-58995955.0. 86 30 1998 *0.849589^6 .0.63877999)  1384 

DATA! 1AI(J)*J=205. 306)=  1185 

*0.62765996, 0,85624997, 0.67961597*0. 79283954, 3, 76052998,0,76896009.  1586 

*0.75699957, 6. 74493958, 0.73266003, 0.72021997, 0.70762998*0. 69453997,  1)97 

*0.6821 6957, 0.66935998 .0.65653993, 0.64374 090 *0.63099998, 0.61 82*995.  1168 

*0.60S42995,S«5«25599950, 57966009, 0,56674999, 0.55385995,0.54101998,  i| 29 

*0,52824998,0.5)555997 ,0.50299996*0. 49052995*  0,478 1 1997. 0.46576995.  1390 

*0.45347995, 0*44 124997, 0.42807995,0. 41696995,0. 40^91995,0.39295996*  31  $1 

*0.38099998, 0,36903995. 0.35*99999.0, 34493495*3. 33291 996, 0,3209939$,  1192 

»0, 30923998 *0-29769999  0. 20643994 *C‘,??5S1 997 ,0.26499999. 0.25*0*997,  1193 

*e.?&49?9«8*0. 23539996*0. 2260$9v9, 0.21699995,0. 20817000,0.18955999,  1394 

*0 .) 81 17999*0.) 8286995*0. 17499895*0, 1672099$, 0*15963995*0, 15228989*  1195 

*0,14515996,0. 13824099*0, 13 hS5997*0«l 2508994* O, 1 1803998  .C.J  3  280*996,  1196 

*0. 10689999. 0.1 0143995*0, 096-12000*  0.Q91C4997.0. 06820000, 0.C815599.S,  119? 

*0.0771 1995*5. 0T285994. 6, 06875956, 0.06481989. 0. *36100000,0. 55735000,  1198 

*O„OS3&SOO0, 0.05060000. 0,04750000, 0,04456000*0. 04178000. 0.03914030,  1199 
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«0. 03663008. 0, 03426000,0.03200000*0. 02990000 ,0.02800000,0.02328000*  3200 

*0. 02473000, 0. 02324000 *0,02188000 *0.02060000 *0.01 938000 rC. 01 81 8000*  2201 

*0.01700000*0.03584700.0,02476800*0.01375606*0,02260800*0.01191800)  1202 

0ATAUA1(J)*J=307.392)=  1203 

*0.01 1084 00*0. 01 030000, 0,00956000 *0.00886200* 0,00820000* 3. 00759000*  1204 

*0*00704600 *0-00655800 *0.006221 00 *0.00572500 *0.00536300 *0.00502800*  1205 

*0. 00472000*0, 00440400*0. 00410000*0. 00361368. 0.00355040.0.00332572*  1206 

8 0. 0031 0320 * 0, 00290938 * 0. 002  .'3080 1 0.0 0256402 *0.0026 0560 *0.00225208*  1207 

*0. 00210000.0.00195415*0.00182068 *0.00 169880.0.00158669. 0.00148312*  1200 

e0*0ei38696.0.031296&6<0, 00121 152*0. 00112966, 0.00105000, 0.00097477,  1209 

«0. 00090656, 0.00084659*0. 00078808* 0*00073625, 0,00068832, 0-00064351*  1210 

®0. 00068104, 0.00056013*0. 00052000, 0.00048184,0. 00044712, 0,00041548,  1221 

•0.00038656,0. 00036000, 0,00033544. 0,00031252*0, 00029080, 0*00027016,  1212 

•0*90*25000.0. 00023102*0. 00021 392, 0, 00019850,0. 00018456. 0, 00017188*  1213 

•0.00016024,0. 00014944*0*  0801 3926, 0.00012954, Os 0001 2000*0.0001 1 103*  1214 

•0,000 13304* 0,00009591 *0.0 0008952 *6. 00008375? 0.00007843* 0.0000 7359*  1 215 

*0.60006896. 0,00006447, 0.00006000,0. 00005572,0. 00005184, 0,00004832*  !2i& 

*0,00004512*0. 00004219,0. 00003948,0. 00003695,0. 00003456, 0,06003226*  1217 

*0,39003000*0.0  )  *218 

END 


r>  no  o  noon 
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SUBROUTINE  EYEBAL  (X.N)  1220 

REV  22MAR71  TO  USE  WITH  BLOCK  DATA  EYES  1221 

PROGRAMMED  BY  8  E  BOUDA.  OCT  1970.  RUN'S  ON  CDC6600.  1222 

A1  IS  EYE  RESPONSE  TABLE.  X  IS  CORRESPONDING  WAVELENGTH.  ANGSTROM  1223 
N  IS  NUMBER  OF  POINTS  .  X<1)=3B00.  XI3921=7710.  1224 

DIMENSION  XS392)  1225 

A1  TABLE  IS  IN  BLOCK  DATA  EYES.  1226 

C0MH0N/B6/ A 1 13921  1227 

N=392  1226 

X <11=3800.  1229 

DO  2  J=2.N  1230 

JH5=J-1  1231 

2  XiJ>=  X( JM1 ) *10,  1232 

VR1TE<6*3>  <X<J1  .Al<Jl*J~l»fil  1233 

3  FORMAT  (/At  6UX.F6.0*E14.7)  11  1234 

NORMALIZE  MAXIMUM  POINT  TO  1.0  1235 

AN=A1<176>  1236 

DO  4  J=1.N  123Y 

4  A1 < J5=AI < J1 /AN  1238 

C  L‘RIT"<6»31  <XU>.A1<J)*J=1«N>  1239 

RET.RN  1240 

END  1241 
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U.S.  Patents 

White  Smoke  Composition ,  3,625,855,  Class  252/305  (7  December  1971). 

Device  for  Producing  White  Smoke  by  Imploding  Red  PJiosphorus , 
3,625,155,  Class  102/90  (7  December  1971),  Coinventor  with  Roscoe 
D-  Dwiggins. 

White  Smoke  Composition  Containing  Red  Phosphorus ,  3,607,472,  Class 
149/19  (21  September  1971). 

Airaxft  Parachute  Slave ,  3,499,385,  Class  102/35  (10  March  1970). 

Illwninating  Flare  Composition  Composed  of  Magnesium 3  Sodium  Nitrate } 
and  a  Vinyl  Terminated  Polysiloxane  Binder ,  3,411,964,  Class  149/19 
(19  November  1963). 

Illuminating  Flare  Composition  Composed  of  Magnesium 3  Sodium  Nitrate , 
and  an  Epoxy  Resin -Polyglycol  Resin  Binder,  3,411,963,  Class  149/19 
(10  November  1968). 

Process  for  Polymerizing  Acrylic  Monomers  with  Strontium  Perchlorate 
for  Pyrotechnics  and  Propellants ,  3,369,946,  Class  149/83  (20  February 
1968). 

Pyrotechnic  Compound  Iris (Glycine)  Strontium  (II)  Perchlorate  and 
Method  for  Making  Same ,  3,296,045,  Class  149/75  (3  January  19£7). 

Plastic  Pyrotechnic  Compositions  Containing  Strontium  Perchlorate 
and  Acrylic  Polymer,  3,258,373,  Class  149/19  (28  June  1966). 


